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ABSTRACT 


A  risk  assessment  was  conducted  to  estimate  the  potential  economic 
losses  through  1993  due  to  the  electrical  effects  of  carbon  fibers 
released  from  U.S.  general  aviation  aircraft.  The  usage  of  carbon  fiber 
(CF)  composites  in  aircraft  structures  is  expected  to  increase  substan¬ 
tially  by  1993.  Aircraft  accidents  often  involve  fires  or  explosions 
that  could  conceivably  release  minute  carbon  fibers,  which  might  disperse 
in  the  atmosphere,  penetrate  buildings  or  enclosures,  and  cause  damaging 
shorts  to  electronic  equipment.  Of  an  estimated  354  annual  general 
aviation  aircraft  accidents  with  fire  in  the  U.S.  in  1993,  approximately 
88  could  involve  aircraft  using  carbon  fibers.  These  accidents  could 
result  in  the  release  of  up  to  1.1  kg.  per  accident,  based  on  forecasts 
of  CF  usage  through  1993  and  experimental  tests  with  burning  CF  composites. 

A  methodology  was  developed  to  compute  estimated  dollar  losses  by 
county  and  equipment  type,  using  a  Poisson  model  for  the  incidence  of 
equipment  failures.  This  approach  incorporated  data  on  the  geographic 
distribution  of  potentially  vulnerable  facilities,  as  well  as  the  mean 
CF  exposure  levels  at  which  various  equipment  would  fail.  The  results 
were  then  statistically  aggregated  to  produce  a  national  risk  profile 
for  estimated  annual  losses  in  1993.  The  expected  national  loss  was 
$253  per  year  (1977  dollars),  and  the  likelihood  of  exceeding  $107,000 
in  annual  losses  was  estimated  to  be  at  most  one  in  ten  thousand.  The 
sensitivity  of  these  results  to  major  input  parameters  was  investigated, 
and  it  was  found  that  under  major  parameter  changes  the  expected  losses 
would  still  remain  low. 


1 


1 .  INTRODUCTION 


1.1  BACKGROUND  AND  OBJECTIVES 

Carbon  fiber  (CF)  composites  are  being  considered  as  an  alternative 
material  in  the  manufacture  of  private  aircraft  because  of  their  light 
weight,  high  strength,  and  design  flexibility.  As  their  production  costs 
decrease,  CF  composites  are  expected  to  find  a  considerable  market  in 
aircraft,  aerospace,  automotive,  and  industrial  applications.  However, 
in  spite  of  the  benefits  of  CF  composites,  a  potential  problem  has  been 
identified  associated  with  the  high  conductivity  of  the  carbon  fibers. 

When  composite  structures  are  exposed  to  fire  of  sufficient  duration  and 
intensity,  it  is  possible  that  the  epoxy  binding  material  will  burn  off, 
releasing  individual  fibers  into  the  atmosphere.  These  fibers,  if 
deposited  on  electronic  equipment,  could  cause  shorts  in  low-voltage 
circuits,  resulting  in  damage  to  the  equipment  and  possible  economic 
losses  for  the  facility  or  community  involved.  The  National  Aeronautics 
and  Space  Administration  (NASA)  has  been  charged  with  the  task  of  inves¬ 
tigating  the  risk  to  the  United  States  as  a  whole  from  potential  releases 
of  CF  in  accidental  aircraft  fires.  As  a  part  of  the  program  of  risk 
assessment,  Arthur  D.  Little,  Inc.,  was  contracted  to  quantify  the  risks 
associated  with  CF  composite  use  in  general  aviation  aircraft  through 
the  year  1993. 

In  order  to  perform  the  risk  assessment,  information  was  gleaned 
from  several  other  agencies  that  are  conducting  parallel  investigations, 
with  NASA  as  the  coordinating  agency.  The  data  incorporated  into  the 
analysis  included  fiber  release  characteristics  for  burning  composites, 
vulnerability  test  results  for  various  categories  of  equipment,  and 
filter  penetration  experiments  which  are  concerned  with  the  ability  of 
single  fibers  to  enter  buildings.  However,  uncertainties  remain  in 
data  inputs  for  certain  crucial  elements  of  the  risk  analysis,  which 
can  introduce  substantial  uncertainty  into  the  magnitude  of  the  resulting 
risk  estimates.  Among  the  areas  of  greatest  uncertainty  are  the  fre¬ 
quency  of  fire  incidents,  the  quantities  of  CF  that  are  actually  released. 
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and  the  equipment-disabling  properties  of  fibers.  In  this  report  we 
have  attempted  to  show  uncertainties  explicitly,  to  make  conservative 
assumptions  where  necessary,  and  to  determine  the  sensitivity  of  our 
risk  estimates  to  these  uncertainties  and  assumptions. 

The  objective  of  the  present  study  was  to  assess  the  1993  national 
risk  of  economic  losses  due  to  the  utilization  of  carbon  fiber  composites 
in  general  aviation  aircraft.  In  formulating  this  objective,  we  identi¬ 
fied  as  sub-objectives  the  projection  of  potential  usage  of  carbon  fiber 
composites  in  U.S.  general  aviation  aircraft  through  1993,  the  develop¬ 
ment  of  an  accident  model  for  general  aviation  aircraft,  the  analysis  of 
of  the  possible  release  amounts  in  general  aviation  accidents,  the  iden¬ 
tification  of  demographic  and  industrial  categories  which  might  be  exposed 
to  such  releases,  and  the  assessment  of  the  economic  consequences  of  a 
given  release.  The  demographic  and  economic  analysis  methods  are  modified 
versions  of  the  approach  used  in  a  parallel  study in  which  we  performed 
a  similar  risk  analysis  for  commercial  air  carriers,  using  Monte  Carlo 
simulation  techniques  to  generate  a  national  risk  profile.  In  the  course 
of  the  present  study,  a  simplified  methodology  was  developed  for  generating 
the  national  risk  profile  by  direct  computation.  This  is  described  in  the 
next  section. 

1.2  METHODOLOGY 

Risk  assessment  of  carbon  fiber  releases  resulting  from  general 
aviation  accidents  is  different  from  previous  risk  assessment  work 
regarding  accidental  CF  releases  from  commercial  aircraft  in  several 
ways.  First,  there  are  substantially  more  general  aviation  accidents 
per  year  than  commercial  aircraft  accidents.  This  difference 
allows  the  utilization  of  analytic  techniques  based  upon  the  statistics 
of  large  numbers.  A  second  difference  is  that  general  aviation  accidents 
are  likely  to  occur  in  flight  or  near  any  of  a  large  number  of  small 
airports,  while  commercial  aircraft  accidents  generally  occur  near  major 
metropolitan  areas.  Finally,  the  most  significant  difference  lies  in 
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the  fact  that  general  aviation  fire  accidents  result  in  relatively  small 
releases  of  carbon  fiber  (compared  to  possible  releases  in  commercial 
aviation),  and  as  a  result  the  failure  probabilities  for  equipment  located 
near  an  accident  are  generally  smaller  than  for  commercial  aviation. 

Because  of  these  differences  an  analytic  model  was  used  for  the 
general  aviation  risk  assessment  instead  of  the  Monte  Carlo  simulation 
procedure  used  for  the  commercial  aviation  analysis.  The  analytic  model 
emphasizes  the  variation  due  to  the  random  nature  of  failure  events, 
rather  than  that  due  to  physical  conditions  such  as  accident  location 
and  weather  variations.  Such  a  model  is  appropriate  when  the  number  of 
failures  per  release  is  low.  In  Appendix  A  it  is  shown  that  since  each 
individual  fiber  or  group  of  fibers  has  a  small  but  finite  probability 
of  causing  a  failure,  and  since  experiments  have  indicated  that  equipment 
failures  obey  an  exponential  probability  law,  then  the  details  of  the 
release  conditions,  with  the  exception  of  the  total  amount  of  fibers 
released,  are  relatively  unimportant.  As  a  result,  given  the  amount 
released,  each  accidental  release  incident  can  be  characterized  by  a 
Poisson  probability  distribution  for  the  number  of  failures.  This 
distribution  can  be  successfully  applied  to  events  for  which  there  are 
a  large  number  of  probabilistic  trials  with  a  low  probability  of  occur¬ 
rence  in  each  trial. 

A  simulation  approach  to  risk  estimation  would  not  be  appropriate 
for  several  reasons.  Since  the  dominant  contribution  to  determining  the 
number  of  failures  is  the  probabilistic  nature  of  the  individual  failures 
(i.e.,  the  Poisson  variation),  the  simulation  approach  requires  a  very 
large  number  of  Monte  Carlo  trials  in  order  to  develop  any  confidence  in 
the  results.  In  addition,  since  general  aviation  accidents  can  occur 
in  widely  dispersed  locations,  a  simulation  would  require  a  data  collec¬ 
tion  effort  that  would  be  prohibitively  costly. 
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The  analytic  methods  that  we  have  developed  for  the  present  applica¬ 
tion  analyze  the  Poisson  failure  process  for  various  equipment  categories 
and  utilize  numerical  calculations  of  probabilities  to  estimate  risk. 

The  analysis  of  equipment  and  facilities  is  performed  on  the  county  level, 
and  the  actual  probability  calculations  are  based  on  mixtures  of  Poisson 
distributions  that  apply  for  each  combination  of  county,  amount  released 
and  equipment  category.  Appendix  A  presents  a  detailed  discussion  of  the 
Poisson  process  approach  and  the  implications  of  low  probability  failures. 

There  were  two  key  parameters  within  the  analytic  model  which 
affected  the  number  of  failures  per  accident.  The  first  was  the  amount 
of  fibers  released  in  an  accident.  By  examining  the  different  types  of 
general  aviation  aircraft  and  their  accident  statistics,  a  distribution 
of  amounts  of  carbon  fibers  potentially  released  in  accidents  was  devel¬ 
oped.  The  second  key  parameter  was  the  density  of  facilities  near  the 
location  of  an  accident.  Thus,  an  important  aspect  of  the  accident 
model  was  a  quantitative  description  of  the  distribution  of  facility 
densities.  The  3,000  counties  in  the  United  States  were  chosen  as  a 
basis  for  estimating  facility  density,  and  hence  a  methodology  was 
developed  to  apportion  accidents  to  the  various  counties. 

The  overall  approach  consisted  of  the  following  steps,  as  illustrated 
in  Figure  1-1 : 

•  A  distribution  of  possible  CF  release  quantities  was 
developed,  based  upon  projected  CF  usage  and  several 
possible  accident  scenarios. 

t  For  each  accident  scenario,  the  surface  integral  of 
exposure  was  estimated.  This  integral  depends  only 
on  the  quantity  of  fibers  released  and  the  fiber 
setting  velocity. 
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FIGURE  1-1  OVERVIEW  OF  METHODOLOGY 
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•  The  conditional  probability  of  a  random  accidental 
fire  occurring  in  each  specific  county  was  estimated 
as  a  function  of  local  and  itinerant  general  avia¬ 
tion  operations. 

•  For  each  county  in  the  U.S.  the  numbers  of  facilities 
in  various  industrial  categories,  as  well  as  private 
residences  and  community  services,  were  enumerated. 
Potentially  vulnerable  equipment  was  identified  within 
each  facility  category. 

•  The  expected  number  of  failures  for  each  class  of 
equipment,  county  location,  and  release  quantity  were 
calculated,  using  information  about  equipment  vulner¬ 
ability  in  terms  of  exposure. 

•  Assuming  that  the  number  of  failures  was  Poisson- 
distributed,  a  probability  distribution  was  generated 
for  the  number  of  failures  per  release  incident, 
aggregated  over  all  counties  and  release  scenarios. 

•  The  proportion  of  failures  occurring  in  each  equipment 
category  was  estimated  and  economic  losses  were  assessed, 
resulting  in  the  statistics  of  dollar  losses  per  release 
incident. 

•  Finally,  the  statistics  of  annual  dollar  losses  were 
obtained  using  the  estimated  total  number  of  fire  incidents 
per  year.  On  the  basis  of  these  statistics,  a  national 
risk  profile  was  generated.  The  national  risk  profile  is 

a  graphical  display  of  the  probability  of  exceeding  various 
levels  of  dollar  loss  as  a  result  of  the  accidental  release 
of  CF  in  a  general  aviation  fire. 
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Chapters  2  to  4  of  this  report  present  the  various  input  data 
required  for  the  risk  analysis,  and  Chapter  5  describes  the  execution 
of  the  above  methodology. 

1,3  RISK  ANALYSIS  PRINCIPLES 

The  concept  of  risk  can  be  defined  as  the  potential  for  realization 
of  unwanted  negative  consequences  of  an  event  or  activity.  In  the  case 
of  this  study,  the  unwanted  negative  consequences  are  the  potential 
economic  losses  due  to  electronic  equipment  failure.  The  event  or 
activity  in  question  is  the  operation  of  general  aviation  aircraft 
utilizing  carbon  fiber  composites.  If  risk  is  due  to  the  presence  of 
some  causative  agent,  such  as  carbon  fibers,  then  the  degree  of  exposure* 
is  measured  by  the  amount  of  that  agent  which  is  potentially  active. 

In  the  past  decade,  an  increasing  amount  of  attention  has  been 
paid  to  problem  areas  involving  activities  with  uncertain  outcomes  which 
might  engender  large  risks.  In  order  to  deal  with  these  problems  the 
field  of  risk  management  has  been  created  and  developed.  Risk  management 
is  a  methodical  scientific  approach  towards  dealing  with  such  risks.  The 
quantitative  aspects  of  risk  management  are  often  referred  to  as  risk 
analysis.  Examples  of  the  application  of  this  approach  are  in  the  areas 
of  nuclear  reactor  safety  and  transportation  of  hazardous  chemicals,  such 
as  liquefied  natural  gases. 

The  practice  of  risk  management  involves  three  basic  steps:  risk 
identification,  risk  measurement,  and  risk  control.  Potential  risks  can 
be  identified  through  experience,  judgment,  or  experimentation.  In  the 
case  of  the  carbon  fiber  problem  the  nature  of  the  risk  is  fairly  well 
understood.  The  major  challenge  lies  in  risk  measurement,  that  is,  in 
determining  the  frequency  of  occurrence  of  events.  Thus,  the  purpose 
of  risk  analysis  is  to  create  an  analytic  framework  permitting  measure- 


In  this  case,  exposure  is  the  time  integral  of  concentration, with  units 
of  fiber-seconds  per  cubic  meter. 
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merit  of  exposure  and  risk.  Finally,  if  the  measured  risk  is  considered 
sufficiently  great,  control  measures  may  be  deemed  necessary.  Control 
measures  would  consist  of  any  modifications  to  the  mechanism  of  risk 
resulting  in  a  reduction  in  the  measured  risk. 

There  are  various  possible  representations  which  can  be  used  to 
quantify  risk.  One  possible  representation  is  the  expected  value  of 
losses  over  a  given  period  of  time.  However,  in  order  to  deal  with 
risks  which  may  fluctuate  over  a  wide  range  of  losses  and  a  corre¬ 
spondingly  wide  range  of  frequencies  of  occurrence,  a  preferred  method 
of  presentation  is  the  risk  profile.  As  discussed  earlier,  a  risk 
profile  is  a  graphical  display  of  risk  showing  the  probability  distri¬ 
bution  for  exceeding  various  levels  of  unwanted  impacts.  A  hypothetical 
example  of  a  risk  profile  is  shown  in  Figure  1-2.  The  activity  in 
question  is  labeled  Activity  1  and  the  risk  profile  for  Activity  1 
shows  that  economic  impact  can  vary  from  $100,000  to  $10  million  with 
probabilities  ranging  from  one  in  a  thousand  to  one  in  ten  thousand. 

This  risk  profile  may  be  compared  against  other  profiles  for  different 
types  of  events,  such  as  the  damage  from  tornadoes.  In  the  diagram, 
two  comparator  risk  profiles  are  shown.  If  risk  control  options  are 
exercised,  it  may  be  possible  to  reduce  the  risk  from  Activity  1  as 
shown  by  the  dotted  curve  at  the  bottom.  The  vertical  lines  are 
confidence  bounds  which  show  the  uncertainty  in  the  estimates  of  risk. 
Even  though  the  actual  risk  may  fall  anywhere  between  these  confidence 
bounds,  the  risk  profile  can  still  be  used  as  an  effective  decision¬ 
making  tool  since  it  both  quantifies  in  an  absolute  sense  the  risks 
imposed  by  Activity  1  and  permits  a  comparison  of  these  risks  relative 
to  other  known  risks. 

1.4  REFERENCES 

1  Arthur  D.  Little,  Inc.,  "An  Assessment  of  the  Risks  Arising  from 
Electrical  Effects  Associated  with  Carbon  Fibers  Released  from 
Commercial  Aviation  Aircraft  Fires".  NASA  Report  No.  NASA  CR-159205 
(February,  1980). 
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2.  CARBON  FIBER  USAGE  FORECASTS 


2.1  INTRODUCTION 

Although  carbon  fiber  (CF)  composite  materials  are  not  presently 
being  used  in  quantity  in  general  aviation  aircraft,  there  is  consider¬ 
able  interest  in  the  potential  utilization  of  these  materials  for  light¬ 
weight,  high  strength  components.^  The  rate  at  which  carbon  fibers 
will  be  introduced  in  the  general  aviation  industry  is  difficult  to 
project  due  to  the  uncertainty  of  prices  and  design  trends.  There  will 
generally  be  modest  amounts  of  carbon  fiber  composites  used  on  general 
aviation  aircraft,  particularly  single  reciprocating  engine  aircraft. 
Eventually  there  may  be  large  amounts  of  carbon  fiber  composites  on  a 
few  private  aircraft.  Lear  Avia,  for  example,  expects  to  be  ultimately 
producing  in  small  quantities  a  jet  aircraft  with  almost  1,000  pounds 
of  carbon  fiber  composites,  although  this  projection  did  not  directly 
affect  our  analysis. 

In  this  section  we  present  an  overview  of  the  projected  usage  of 
carbon  fiber  composites  in  general  aviation,  and  then  estimate  the 
amount  of  CF  in  the  1993  fleet.  These  projections  are  based  mainly  on 
discussions  between  Arthur  D.  Little,  Inc.,  aircraft  manufacturers,  and 
NASA,  and  represent  only  an  approximate  forecast. 

2.2  OVERVIEW  OF  COMPOSITE  USAGE 

The  general  aviation  and  business  aircraft  industry,  which  includes 
both  fixed  wing  and  rotary  aircraft,  continues  to  represent  a  very  small 
portion  of  fiber  consumption  within  the  aerospace  industry.  In  1978 
this  sector  consumed  no  more  than  3,000  kg.  of  carbon  fibers  which  repre¬ 
sented  less  than  4%  of  total  aerospace  usage.  The  major  portion  of  this 
volume,  probably  as  much  as  65%  to  70%,  was  used  by  rotary  aircraft 


^Larson,  G.C.,  Composite  Materials  and  General  Aviation,  Business  and 
Commercial  Aviation  (September,  1979). 
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manufacturers,  with  the  remainder  split  among  eight  North  American  fixed 
wing  aircraft  producers. 

The  rotary  aircraft  (helicopter)  industry  is  made  up  of  the  follow¬ 
ing  five  companies:  Sikorsky,  Bell -Helicopter  Textron,  Boeing-Vertol , 
Hughes, and  Kaman.  The  major  output  of  these  companies,  currently  esti¬ 
mated  at  65%  to  70%,  goes  into  military  service.  The  remainder  is  sold 
to  various  other  customers  in  the  commercial  and  general  aviation  market. 
All  five  companies  are  currently  using  carbon  fiber  to  varying  degrees, 
both  in  research,  development,  and  evaluation  programs  sponsored 
by  NASA  or  other  government  groups,  and  in  selected  operational  appli¬ 
cations. 

The  helicopter  producers  have  generally  made  only  limited  inroads 
in  the  use  of  CF  epoxy  composites  relative  to  the  commercial  and 
military  fixed  wing  aircraft  sectors,  and  in  most  cases  on  a  retrofit 
basis.  The  primary  reasons  for  this  have  been  the  limited  number  of 
units  produced  over  which  the  development  and  production  costs  can  be 
spread,  and  the  general  acceptability  of  considerably  less  expensive 
materials  in  most  applications.  It  appears  the  contamination  issue 
has  not  been  a  deterring  factor  in  the  use  of  CF.  The  primary  appli¬ 
cation  for  CF  has  been  as  a  stiffening  agent  in  such  areas  as  tail  rotor 
spars,  spar  blades,  main  rotors,  fairings,  and  horizontal  stabilizers. 

In  most  cases  where  composite  materials  are  utilized,  fiber  glass  and/or 
Kevlar  constitute  the  major  portion  of  the  fiber  requirement  with  CF 
used  in  limited  quantities  for  local  stiffening.  Normal  ratios  for 
glass  and/or  Kevlar  versus  graphite  range  between  5:1  to  10:1. 

On  a  weight  basis,  total  carbon  fiber  consumption  within  the  heli¬ 
copter  industry  in  1978  approached  7,000  kg.  Only  2,000  of  these  kg. 
were  believed  to  be  used  in  areas  other  than  military  and  commercial 
applications,  with  the  military  sector  accounting  for  the  lion's  share. 

On  a  company-by-company  basis,  Sikorsky  represented  the  largest  consumer, 
utilizing  about  4,000  of  the  total  7,000  kg.  The  vast  majority  of  this 
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material,  however,  went  into  military  aircraft  with  the  S-76  "Spirit" 
being  the  sole  general  aviation  aircraft  to  utilize  CF.  The  next 
largest  consumers  were  Boeing-Vertol  and  Kaman,  each  utilizing  about 
800  kg.  of  fiber.  Here  again  the  major  portion  went  toward  military 
applications.  The  remaining  industry  usage  was  split  between  Bell,  Hughes 
and  various  subcontractors  for  the  industry,  such  as  Fiber  Science. 

The  predominant  application  for  CF  to  date  has  been  in  rotor  and 
spar  blades  for  the  retrofit  market.  Uni-directional  tape  has  been  the 
predominant  material  form  used.  Some  woven  material  is  going  into  spar 
booms,  but  other  than  this,  woven  product  usage  is  minimal.  Current 
evaluation  programs  are  under  way  at  the  various  manufacturers  for  com¬ 
ponents  such  as  rotor  hubs,  transmission  shafts,  gearbox  struts,  landing 
gear  struts  and  various  housings.  It  appears  likely  that  uni -directional 
tape  will  be  the  predominant  material  used  for  these  applications  as  well; 
however,  woven  material  is  expected  to  be  used  to  at  least  some  degree. 

In  terms  of  forecast  carbon  fiber  consumption  in  non-military  heli¬ 
copter  applications,  we  estimate  usage  will  experience  a  steady  climb  up 
to  a  level  of  12,000  to  15,000  kg.  by  1985,  which  represents  at  least  a 
six-fold  increase  over  current  levels.  Probably  as  much  as  30%  of  this 
usage  will  go  into  research  and  development  and  new  component  evaluation 
programs,  with  the  remainder  heavily  weighted  toward  current  applications. 
Uni-directional  tape  will  continue  to  comprise  the  lion's  share  of  the 
product  form,  with  woven  material  gaining  only  slightly  in  popularity. 

It  is  difficult  to  predict  usage  levels  after  this  period  since  any 
number  of  applications  currently  under  evaluation  might  be  specified 
during  this  period,  any  of  which  could  greatly  affect  the  total  consump¬ 
tion  levels.  We  feel,  however,  that  fiber  glass  and  Kevlar  will  continue 
to  represent  the  major  portion  of  composite  material  consumed  and  CF  will 
continue  to  be  used  only  in  local  stiffening  applications  in  conjunction 
with  the  above  materials. 

With  respect  to  the  business  and  general  aviation  fixed  wing  air¬ 
craft  industry,  the  primary  manufacturers  include  Lear  Avia,  Beech,  Gates 
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Lear jet,  Cessna,  Grumman,  Piper,  and  Rockwell  International  (all  of  the 
United  States)  and  Canadair  and  deHavilland  of  Canada.  Only  a  few 
of  the  above  companies,  namely  Lear  Avia,  Beech,  Cessna  and  Gates  Learjet, 
are  currently  using  carbon  fiber,  predominantly  in  evaluation  programs. 

In  addition  to  these  aircraft  manufacturers,  there  are  a  number  of  other 
small  companies  using  limited  amounts  of  carbon  fiber  in  components  made 
for  both  domestic  and  foreign  aircraft  producers. 

As  was  the  case  in  the  helicopter  industry,  the  use  of  CF  has  been 
considerably  limited  by  low  production  volumes  and  the  acceptability  of 
more  conventional  and  considerably  less  expensive  materials.  CF  has 
usually  only  been  specified  as  a  local  stiffening  agent  for  components 
produced  predominantly  from  fiber  glass  at  Kevlar,  the  two  largest  of 
which  are  spoilers  and  flaps.  Components  produced  entirely  from  fiber 
glass  and/or  Kevlar  include  such  products  as  wing  tips,  wing/fuselage 
spars,  radomes,  engine  cowlings,  tail  covers  and  various  fairings. 

The  one  major  exception  to  the  limited  appl icational  usage  for  CF 
is  Lear  Avia,  where  they  have  been  developing  an  all  composite  aircraft. 

This  aircraft,  labeled  the  Learfan,  will  be  made  entirely  of  carbon/epoxy 
and  other  composite  materials  with  the  exclusion  of  the  engine,  windows, 
and  assorted  other  subcomponents.  The  total  empty  weight  of  this  seven 
to  eight  passenger,  turboprop  aircraft  will  be  1,650  kg.,  700  kg.  of  which 
will  be  comprised  of  composite  materials. 

The  Learfan  will  use  both  pure  fiber  laminates  as  well  as  hybridized 
forms,  including  CF,  glass  and  Kevlar.  Possibly  even  boron  will  be  used 
as  well  as  phenol ics  and  polymides  where  higher  temperatures  and/or  stresses 
present  a  problem.  Woven  fabric  is  expected  to  comprise  the  majority  of 
the  material,  although  liberal  amounts  of  uni -directional  and  multi -direc¬ 
tional  tape  will  also  be  employed. 

This  design  is  a  radical  deviation  from  the  conventional  applications 
of  CF  in  strictly  secondary  structural  components.  Lear  hopes  to  have  this 
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aircraft  test  flown  within  one  year  and  will  then  go  into  full  production 
shortly  afterwards,  expecting  to  attain  a  production  rate  of  250  to  300 
per  year  by  1983.  While  there  is  still  a  great  deal  of  market  speculation 
about  its  flight  capabilities,  potential  CF  contamination  issues,  and 
selling  price  (supposedly  the  Learfan  will  be  priced  at  about  $1.2  million), 
Lear  states  it  already  has  orders  for  over  110  aircraft.  The  conductive 
filament  problem  has  not  been  resolved  to  date  and  this  might  affect  FAA 
certification.  Because  of  the  uncertainty  concerning  the  Learfan  and  the 
limited  amount  of  total  composite  weight  involved,  compared  to  the  1993 
forecast,  we  did  not  consider  its  impact  in  the  analysis. 


Of  the  remaining  business  and  general  aviation  aircraft  producers, 
only  Cessna  and  Beech  appear  to  be  using  CF  and  only  in  extremely  limited 
amounts.  Cessna,  with  almost  50%  of  the  general  aviation  aircraft  market, 
is  the  only  one  of  the  two  currently  using  CF  in  an  operational  applica¬ 
tion.  The  material  is  used  as  a  stiffening  agent  in  the  spoilers  and 
flaps  of  the  new  Citation  III  aircraft.  Beech  is  still  evaluating 
various  components  and  to  date  has  elected  not  to  specify  CF. 

The  current  consumption  level  of  carbon  fiber  by  Lear  Avia  and  the 
others  was  believed  to  be  no  higher  than  1,000  kg.  in  1978.  Lear  Avia 
alone  accounted  for  well  over  700  kg.  of  this.  While  it  is  uncertain  at 
this  time  what  kind  of  carbon  fiber  consumption  levels  will  be  attained 
at  Lear  Avia,  depending  upon  the  success  of  failure  of  the  Learfan,  it 
is  highly  unlikely  that  general  consumption  will  increase  very  rapidly 
over  the  next  five-year  period.  The  primary  reason  for  this  pessimistic 
forecast  is  the  high  cost  premiums  associated  with  CF  and  the  questionable 
need  for  such  a  high  performance  material  in  general  aviation  application. 
Even  with  a  growth  rate  similar  to  that  in  rotary  aircraft  applications, 
carbon  fiber  usage  within  this  sector,  excluding  Lear  Avia,  will  be  no 
larger  than  1,500  to  2,000  kg.  by  1985,  which  is  insignificant  relative 
to  total  consumption. 

Lear  Avia  could  greatly  affect  the  consumption  levels  of  this  sec¬ 
tor;  however  even  with  the  most  optimistic  growth  rates  possible,  total 
usage  would  still  be  small  relative  to  commercial  and  military  usage. 
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Lear  Avia  projects  that  by  1985  annual  graphite  fiber  requirements  could 
be  in  excess  of  30,000  kg. 

The  overall  future  of  graphite  fiber  within  the  general  aviation 
sector  therefore  is  not  expected  to  be  very  promising.  Current  usage 
is  limited,  and  even  with  significant  increases  in  applications,  because 
of  the  small  production  volume  relative  to  the  commercial  and  military 
sectors,  this  area  will  continue  to  account  for  less  than  5%  of  total 
aerospace  fiber  consumption.  We,  nevertheless,  have  used  conservatively 
high  projections  of  1993  carbon  fiber  usage,  as  described  in  Section 
2.4  below. 

2.3  GROWTH  RATES  FOR  GENERAL  AVIATION  AIRCRAFT 

For  the  purposes  of  this  study,  we  have  classified  general  aviation 
aircraft  into  three  major  categories: 

1.  Single  reciprocating  engine  craft 

2.  Multi-engine  and  jet  aircraft,  and 

3.  Helicopters,  non-fixed  wing  and  non-powered  aircraft 

The  identification  of  these  three  categories  was  based  upon  the  accident 
analysis  presented  in  Chapter  3,  and  upon  size  and  structural  differences 
between  the  different  types  of  general  aviation  aircraft.  Of  the  third 
category,  helicopters  represent  the  overwhelming  majority  of  aircraft. 

The  first  two  categories  represent  fixed  wing  powered  aircraft.  Of 
these  type  of  aircraft,  the  first  category  represents  larger  aircraft 
and  is  expected  to  use  larger  amounts  of  carbon  fiber  composites  in  the 
future. 

Historical  growth  rates  for  the  different  classes  of  aircraft  are 
presented  in  Table  2-1.  As  noted,  the  growth  rates  for  the  three 
classes  of  aircraft  are  4.8%,  6.3%,  and  8.6%  per  year  respectively. 

Since  the  three  classes  of  aircraft  are  used  to  define  accident 
categories,  the  three  growth  rates  were  applied  in  Chapter  3  to  estimate 
the  1993  conditional  probability  that  an  accident  involves  an  aircraft 
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TABLE  2-1 


ACTIVE* GENERAL  AVIATION  AIRCRAFT 


Fixed  Wing  Aircraft 


Year 

Multi-Engine 

Single-Engine 

Rotor 

1976 

25,684 

144,941 

6408 

1975 

24,559 

136,651 

5200 

1974 

23,418 

131 ,932 

4224 

1973 

21,929 

126,217 

4213 

1972 

19,849 

120,446 

4100 

1971 

17,855 

109,256 

3916 

Growth  Rate  (%/yr) 

6.3% 

4.8% 

8.6% 

*active  aircraft  is  defined  as  one  which  has  a  current  registration 
and  was  flown  during  the  previous  calendar  year 


Source:  Aerospace  Facts  and  Figures,  1978/79 
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of  a  given  type.  These  growth  rates  also  imply  that  the  number  of  air¬ 
craft  in  use  in  1993  will  be  321,000,  72,500  and  25,000,  respectively. 

2.4  PROJECTED  CF  USAGE  IN  1993 

Any  future  projections  of  carbon  fiber  usage  in  general  aviation 
aircraft  are  subject  to  some  uncertainty.  However,  based  on  discussions 
with  manufacturers  it  is  anticipated  that  approximately  25%  of  the 
general  aviation  fleet  will  be  using  carbon  fibers  by  1993.  This  pro¬ 
jection  is  based  on  the  large  increases  in  fleet  size  between  now  and 
1993.  A  significant  number  of  the  newly  constructed  planes  will  be 
built  with  carbon  fiber  composites.  Because  of  the  uncertainty  involved 
in  these  projections,  our  risk  profiles  were  subjected  to  sensitivity 
analyses  (see  Chapter  6). 

Based  on  these  projections  there  will  be  80,250  single  reciprocating 
engine  craft,  18,125  multi-engine  and  jet  aircraft,  and  6,250  helicopter 
and  non-powered  aircraft  using  carbon  fibers  in  1993.  It  is  anticipated 
that  there  will  be  not  more  than  1.3  million  kilograms  of  carbon  fiber 
composites  manufactured  for  use  on  general  aviation  aircraft  by  1993. 

Using  this  conservative  forecast,  and  based  on  the  belief  that  there  will 
be  large  amounts  used  on  jet  aircraft  and  still  larger  amounts  on  heli¬ 
copters,  we  used  the  following  projections  for  amount  of  composite  per 
aircraft  carrying  CF  in  1993. 

Single  reciprocating  engines  7  kilograms 

Multi-engine  and  Jet  Aircraft  20.5  kilograms 

Helicopter,  non-fixed  wing 
and  non-powered  50.5  kilograms 

These  projections  were  used  in  Chapter  3  to  determine  the  possible  release 
amounts  in  general  aviation  accidents. 
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3, _ PROJECTION  OF  ACCIDENTAL  CARBON  FIBER  RELEASES 


3.1  INTRODUCTION 

An  important  aspect  of  the  risk  assessment  was  a  projection  of  the 
frequency  and  type  of  general  aviation  accidents  that  would  result  in 
accidental  releases  of  carbon  fibers.  This  involved  an  analysis  of 
general  aviation  accident  histories  from  NTSB  data  and  characterization 
of  the  conditions  associated  with  potential  CF  releases.  Due  to  the 
complexity  of  the  historical  accident  data,  the  analytical  details  are 
presented  in  Appendix  B.  This  chapter  discusses  only  the  highlights  of 
the  analysis  and  presents  the  results  that  were  subsequently  used  for 
risk  estimation. 

3.2  GENERAL  AVIATION  ACCIDENT  CHARACTERISTICS 

Three  categories  of  general  aviation  aircraft  were  defined,  as 
mentioned  in  Chapter  2.  It  was  determined  that  accidents  for  these 
different  aircraft  types  could  be  characterized  to  a  sufficient  extent 
by  two  key  variables  -  the  phase  of  operation  and  the  level  of  damage 
sustained.  Due  to  the  Poisson  methodology  used  in  this  risk  assessment, 
weather  conditions  and  precise  accident  locations  were  not  significant 
in  estimating  the  number  of  equipment  failures  per  CF  release  incident. 

The  phase  of  operation  was  classified  as  either  cruise  or  on  or 
near  airport,  the  latter  classification  encompassing  takeoffs,  landings 
and  static  or  taxi  phases.  The  level  of  damage  was  classified  as  either 
substantial  damage  or  total  destruction,  in  accordance  with  NTSB  designa 
tions.  The  conditional  probability  of  having  a  specific  aircraft  type, 
phase,  and  damage  level  are  shown  in  Table  3-1.  These  probabilities 
represent  the  fraction  of  incidents  in  which  each  combination  of  charac¬ 
teristics  would  occur.  For  example,  given  that  an  accident  occurs, 
there  is  a  one-third  chance  that  it  will  involve  total  destruction  of 
a  single  reciprocating  engine  craft  during  cruise.  The  figures  in 
Table  3-1  are  based  on  historical  data  and  are  adjusted  to  reflect  the 
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TABLE  3-1 

CONDITIONAL  PROBABILITIES  OF  ACCIDENT  CHARACTERISTICS  IN  1993 


Cruise -  — On  or  Near  Airport — 


Total 

Destruction 

Substantial 

Damage 

Total 

Destruction 

Substantial 

Damage 

Non-Fixed  or 
Non-Powered 

.072 

.013 

.043 

.014 

Single  Reciprocating 

.333 

.023 

.203 

.034 

Multiple  or  Jet 

.101 

.014 

.122 

.028 
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growth  rates  for  each  aircraft  type  through  1993,  which  were  presented 
in  Chapter  2. 

To  perform  the  national  risk  calculation,  it  was  necessary  to 
estimate  the  total  number  of  general  aviation  accidents  in  the  U.S.  in 
1993,  and  to  allocate  these  accidents  to  different  counties  according 
to  the  degree  of  air  traffic  activity.  The  historical  accident  rate 
appeared  to  be  relatively  constant,  and  hence  was  projected  to  remain 
stable  through  1993,  Consequently,  the  annual  number  of  general  aviation 
accidents,  including  commuter  and  air  taxi  operations,  was  estimated  at 
354,  Of  these  about  25%,  or  88  accidents,  would  involve  CF  according  to 
the  usage  projections  in  Chapter  2.  These  accidents  were  allocated  to 
various  counties  based  on  the  estimated  number  of  local  and  itinerant 
operations  in  each  county.  The  details  of  this  procedure  are  presented 
in  Appendix  B. 

3.3  PROJECTION  OF  AMOUNT  RELEASED 

Given  that  an  accident  occurs  involving  a  general  aviation  aircraft 
using  carbon  fiber  composites  and  that  a  subsequent  fire  results,  a  po¬ 
tential  exists  for  release  of  carbon  fibers.  In  order  to  estimate  the 
resulting  damage,  it  is  necessary  to  know  the  potential  exposure  to  the 
surrounding  area  due  to  carbon  fibers.  The  phenomenon  of  carbon  fiber 
release  and  dispersion  involves  a  complex  chain  of  events,  and  to  physi¬ 
cally  model  these  events  would  require  a  knowledge  of  the  fire  parameters 
such  as  pool  size,  duration,  and  amount  of  fuel  burned,  as  well  as  the 
weather  conditions  at  the  time  of  accident.  Since  these  parameters  would 
be  difficult  to  specify  in  the  case  of  randomly  located  general  aviation 
accidents,  we  have  adopted  a  simplified  methodology  (as  described  quali¬ 
tatively  in  Chapter  1  and  in  detail  in  Appendix  A)  which  circumvents  the 
need  for  most  information.  The  only  information  necessary  is  the  total 
amount  of  carbon  fibers  released  in  the  fire,  since  the  distribution 
of  the  number  of  failed  equipment  becomes  independent  of  the  other  release 
conditions.  This  section  presents  the  assumptions  concerning  the  amount 
of  carbon  fibers  released  in  a  general  aviation  fire  accident,  which 
depends  on  two  factors: 
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•  The  amount  of  CF  on  the  aircraft 

•  The  severity  of  the  fire  and/or  explosion 

It  was  assumed  that  the  amount  released  will  be  proportional  to  the 
total  amount  of  carbon  fiber  composites  being  used  on  the  aircraft.  As 
noted  in  Chapter  2,  there  are  three  types  of  general  aviation  aircraft, 
namely  single  reciprocating  engine  planes,  multiple  and  jet  engine  planes, 
and  non-fixed  wing  and  non-powered  planes.  The  amounts  of  carbon  fibers 
being  used  on  these  craft  are  presented  in  Chapter  2.  It  is  forecasted 
that  25%  of  the  1993  fleet  will  be  using  carbon  fibers  and  hence  the 
probability  that  an  accident  will  result  in  a  CF  release  is  25%. 

The  model  for  general  aviation  aircraft  accidents  described  in 
Section  3.2  defines  severity  in  terms  of  the  NTSB  classification  of 
substantial  damage  and  total  destruction.  These  classifications  were 
used  for  the  purposes  of  determining  the  percentage  of  fibers  released 
in  an  accident.  Based  on  experimental  findings  reported  in  the  analysis 
of  commercial  aviation  aircraft*  it  was  estimated  that  not  more  than 
1%  of  carbon  fibers  would  be  released  in  a  fire  and  that  not  more  than 
2.5%  would  be  released  in  most  fire  and  explosion  scenarios.  In  addition, 
in  commercial  aviation  accidents  not  all  of  the  carbon  fiber  composites 
would  be  consumed  in  a  fire,  while  the  explosive  mode  represents  only  a 
small  minority  of  commercial  aviation  accidents. 

Based  on  these  considerations,  we  conservatively  assumed  that  the 
carbon  fiber  released  in  a  general  aviation  accident  would  be  2%  and 
0.5%  of  the  total  carbon  fiber  composite,  respectively,  for  total 
destruction  and  substantial  damage  accidents.  Since  the  carbon  fiber 
mass  comprises  approximately  70%  of  the  mass  of  the  carbon  fiber  compo¬ 
site,  these  assumptions  represent  2.9  and  0.7%  of  the  carbon  fiber  mass 
respectively  for  total  destruction  and  substantial  damage  accidents. 


*See  Reference  1 ,  Chapter  1 


22 


From  the  assumptions  in  Chapter  2  about  composite  utilization  in  1993, 
we  determined  the  total  mass  of  carbon  fibers  released  for  each  combina¬ 
tion  of  accident  severity  and  aircraft  type.  These  data  are  presented 
in  Table  3-2.  The  maximum  possible  release  is  1.09  kilograms  for  the 
case  of  a  total  destruction  helicopter  accident. 
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TABLE  3-2 


MASS  OF  CARBON  FIBERS  RELEASED  (KG.)  FOR  GENERAL 
AVIATION  ACCIDENTS  -  1993 


Single  Reciprocating 
Engi ne 


Multiple  or  Jet  Engine 


Non-Fixed-Wing  or 
Non-Powered 


Substantial 

Damage 

Acci dent 

Total 

Destruction 

Accident 

.034 

.14 

.10 

.41 

.27 

1.09 
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4.  DEMOGRAPHIC  ANALYSIS  OF  VULNERABLE  FACILITIES 


4.1  INTRODUCTION 

The  national  risk  profile  for  economic  losses  resulting  from  acci¬ 
dental  carbon  fiber  releases  from  general  aviation  aircraft  was  based 
on  the  distribution  of  facilities  with  vulnerable  equipment.  A  set  of 
parameters  was  selected  to  describe  each  U.S.  county  for  the  purposes 
of  the  risk  analysis  presented  in  Chapter  5.  These  parameters  pertain 
to  demographic  data  which  are  readily  available  from  published  sources. 
This  chapter  presents  the  basis  for  the  demographic  analysis,  as  well 
as  the  economic  analysis  of  the  consequences  of  failures.  Most  of  the 
data  utilized  here  were  developed  in  a  parallel  study  of  air  carrier 
fires  and  CF  releases.  (See  Reference  1,  Chapter  1). 

4.2  METHODOLOGY 

The  first  step  in  the  analysis  was  to  represent  the  facilities 
considered  to  be  potentially  vulnerable  by  demographic  classes  such 
as  households  or  the  Standard  Industrial  Classification  (SIC)  codes  for 
businesses.  For  several  other  facility  categories,  indices  were 
required  where  actual  data  on  facilities  were  not  available;  for 
instance,  population  was  used  as  a  surrogate  to  measure  the  amount  of 
police  and  fire  protection  services.  Table  4-1  shows  the  facility 
categories  and  the  demographic  data  used  to  represent  each  facility 
category.  Table  4-2  shows  the  data  sources  for  each  of  the  demographic 
data  classifications.  The  assignment  of  facility  categories  studied  in 
the  economic  analysis  to  demographic  data  categories  involved  some 
aggregation.  For  example,  the  general  manufacturing  category  includes 
equipment  classes  identified  in  specific  manufacturing  environments 
which  were  taken  as  representative  of  the  level  of  vulnerable  equipment 
in  all  manufacturing  plants. 

Given  the  data  categories  for  facilities,  the  amount  of  activity 
in  terms  of  number  of  pieces  of  equipment  in  each  county  was  determined 
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TABLE  4-1 


FACILITY  AND  DEMOGRAPHIC  CATEGORIES 


Households 

Police  Protection  Services 
Fire  Protection  Services 
Post  Office  Sorting  Centers 
Subways 

Commuter  and  Intercity 
Rail  road 

General  Manufacturing 

Manufacturers  of  Electronic 
Equipment 

Telephone  Company  Switching 
Facilities 

Radio  and  Television 
Broadcasting 

General  Merchandise  Retailers 

Retail  Grocers 

Financial  and  Insurance 
Services 

Computer  Services 

Electronic  R&D  Firms  and 
Universities 

Hospitals 

Airport  Services 

Automobile  and  Truck 
Assembly 


Families 

Population 

Population 

Population 

Number  of  Rapid  Transit  Vehicles 


Railroad  Terminals 
SIC  Code  19 


SIC  Codes  3573,  3650,  3660,  3670 
Families 


SIC  Codes  4830,  4890 

SIC  Codes  5310,  5600,  5700,  5900 

SIC  Code  5410 


SIC  Codes  6020,  6100,  6200,  6300 
SIC  Code  7370 


SIC  Codes  7391,  8220 

Number  of  Hospital  Beds 

Number  of  Air  Carrier  Operations  -  1977 

SIC  Code  3710 
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TABLE  4-2 


DEMOGRAPHIC  DATA  SOURCES 


Demographic  Data  Category  Data  Source 


SIC  Data 


U.S.  Census  Bureau,  1976  County  Business 
Patterns 


Families,  Population,  U.S.  Census  Bureau,  1977  County  and  City 

Number  of  Hospital  Beds  Data  Book 


Number  of  Rapid  Transit  American  Public  Transit  Association 

Vehicles 


Railroad  Terminals  The  Official  Railway  Guide,  North 

American  Passenger  Travel  Edition, 
July/August  1979 


Number  of  Air  Carrier  U.S.  Department  of  Transportation, 

Operations  -  1977  Federal  Aviation  Administration, 

Terminal  Area  Forecasts,  Fiscal  Years 
1979-1990 
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from  scaling  factors.  These  scaling  factors  included  such  measures  as  the 
average  number  of  employees  per  SIC  category.  For  each  facility  surveyed 
in  the  economic  analysis,  the  number  of  pieces  of  equipment  and  the  value 
of  the  scaling  factor  for  that  facility  were  determined.  From  the  survey, 
standard  factors  were  developed,  such  as  one  piece  of  equipment  class  x  for 
every  1,000  employees  in  SIC  category  y.  In  this  manner,  the  number  of 
pieces  of  equipment  in  each  category  of  vulnerable  equipment  was  determined 
for  each  facility  category.  Appendices  C  and  D  contain  listings  of  the 
equipment  categories,  with  the  scaling  factors  used  for  each  facility 
type. 


For  every  facility  and  equipment  combination,  the  following  parameters 
were  identified:  the  mean  dosage  for  failure,  the  transfer  functions  for 
outside  to  inside  CF  exposure,  and  the  dollar  cost  per  failure.  For  con¬ 
venience  in  the  risk  computation,  described  in  detail  in  Chapter  5,  the  mean 
dosage  for  failure  and  the  transfer  functions  were  combined  to  develop  the 
effective  mean  outside  dosage  F  for  failure.  When  there  was  a  range  of 
transfer  functions  depending  on  building  characteristics,  the  arithmetic 
mean  of  the  high  and  low  transfer  functions  was  used;  this  procedure  resulted 
in  a  number  of  about  the  same  order  of  magnitude  as  the  high  end  of  the 
transfer  function  range,  which  is  a  consistently  conservative  assumption. 
Equipment  categories  which  had  equivalent  E  values  and  equivalent  demo¬ 
graphic  data  categories  were  combined  for  efficiency  in  computer  process¬ 
ing.  The  dollar  cost  per  failure  of  one  piece  of  equipment  was  derived 
as  the  weighted  average  of  the  unit  costs  for  each  equipment  category. 

Given  the  estimate  of  the  number  of  pieces  of  equipment  for  each 
facility  category  and  equipment  type,  the  computer  procedure  described 
in  Chapter  5  could  be  implemented,  providing  probabilities  of  equipment 
failure  for  each  category.  The  risk  profile  for  dollar  losses  was 
derived  by  combining  these  probabilities  with  the  dollar  loss  per  failure 
of  equipment.  These  losses  were  taken  as  the  sum  of  the  equipment  repair 
and  facility  disruption  costs  per  failure  of  equipment.  In  theory,  this 
procedure  could  overestimate  losses  if  the  expected  number  of  pieces  of 
equipment  failing  in  a  single  facility  were  greater  than  one;  in  that 
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case  the  facility  disruption  cost,  which  might  not  increase  beyond  the 
first  equipment  failure,  would  be  overestimated.  However,  with  the  CF 
releases  being  very  low  relative  to  the  E  values,  the  expected  number 
of  equipment  failures  in  any  facility  would  always  be  systematically  lower 
than  one,  AppendixC  shows  the  estimated  dollar  losses  per  equipment  failure. 

4.3  SUMMARY  OF  ECONOMIC  ANALYSIS 

Details  of  the  economic  loss  analysis  may  be  found  in  Reference  1, 
Chapter  1.  Some  of  the  major  observations  that  resulted  were  as  follows: 

•  Most  industrial  and  commercial  facilities  are  equipped 
for  repair  or  replacement  of  electronic  devices  in  the 
event  of  failures  during  normal  operation. 

•  Equipment  which  are  critical  to  the  operation  of  a 
facility,  such  as  computers,  are  usually  given  special 
protection,  and  backup  procedures  are  often  available 
to  prevent  facility  shut-down, 

•  Much  of  the  electronic  equipment  examined  is  virtually 
invulnerable  to  the  expected  levels  of  indoor.  CF  ex¬ 
posure  (at  most  10^  fiber-seconds/m. 3)  due  to  protective 
cabinets  and  filtration  systems. 

Consequently,  there  are  few  instances  in  which  a  facility  would 
experience  significant  economic  losses  as  a  result  of  CF  exposure.  The 
maximum  dollar  loss  estimated  for  a  single  equipment  failure  was  the 
$65,800  attributed  to  the  loss  of  a  transformer  substation  switch.  The 
generally  low  failure  costs  are  reflected  in  the  low  risk  estimates 
derived  in  Chapter  5. 
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5.  DEVELOPMENT  OF  NATIONAL  RISK  PROFILE 


5.1  INTRODUCTION 

This  chapter  describes  the  methodology  used  to  determine  the 
national  risk  profile  and  presents  an  interpretation  of  the  results. 

The  methodology  utilized  a  computer  model  to  calculate  the  probability 
distribution  for  the  consequences  of  a  single  accident.  These  single 
accident  results  were  then  extrapolated  to  obtain  a  national  estimate 
of  expected  annual  losses. 

The  remainder  of  the  chapter  is  divided  into  four  sections. 

Section  5.2  presents  the  methodology  and  results  for  the  potential 
economic  losses  in  a  single  automobile  accident.  The  mathematical 
basis  for  the  methodology  in  this  section  is  presented  in  Appendices 
A  and  E.  In  Section  5.3,  the  results  for  a  single  incident  are 
extrapolated  to  an  annual  risk  profile.  The  extrapolation  technique 
uses  the  distribution  of  dollar  losses  in  a  single  accident  to  derive 
an  annual  dollar  loss  distribution  based  on  an  expected  88  accidents 
per  year  involving  CF  composites.  In  Section  5.4,  results  of  a  sensi¬ 
tivity  analysis  are  presented.  It  is  noted  that  the  change  in  annual 
dollar  loss  probabilities  with  respect  to  the  changes  in  input  parameters, 
such  as  release  amounts,  can  be  represented  by  a  very  simple  mathematical 
relationship.  Finally  Section  5.5  contains  a  summary  discussion  of  the 
results. 

Before  examining  the  details  of  the  methodology  it  is  important 
to  understand  the  principle  of  the  Poisson  approach.  For  a  given  release 
scenario  and  equipment  type,  the  number  of  equipment  failures  may  be 
approximated  by  a  Poisson  distribution.  The  mean  number  of  failures  is 
given  by  integrating  the  equipment  density  over  the  area  in  question  and 
multiplying  by  the  equipment  failure  probability,  which  is  nearly  linear 
in  E  for  low  values  of  the  exposure  E.  Under  modest  assumptions,  we  can 
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aggregate  over  many  release  scenarios,  and  show  that  the  average  number 
of  failures  is  proportional  to  the  surface  integral  S  of  the  exposure, 
which  in  turn  may  be  shown  to  depend  only  on  the  amount  released  and  the 
fiber  settling  velocity.  Thus  an  expression  is  obtained  for  the  mean 
number  of  failures  per  incident  in  terms  of  just  the  average  facility 
density,  the  amount  of  CF  released,  and  the  equipment  vulnerability. 

This  enables  us  to  assess  the  risk  without  requiring  detailed  data  on 
accident  conditions  or  geographic  locations. 

Most  of  the  technical  details  of  the  methodology  are  presented  in 
the  appendices.  There  are,  however,  some  fundamental  mathematical 
relationships  that  control  the  results  developed  in  this  report.  These 
relationships  are  presented  below  to  emphasize  their  importance  in  the 
final  analysis.  A  glossary  of  symbols  used  in  the  relationships  dis¬ 
cussed  in  this  chapter  is  presented  in  Table  5-1. 

The  first  key  relationship  is  between  X,  the  expected  number  of 
equipment  failures  in  an  accident,  and  such  parameters  as  the  amount  of 
carbon  fibers  released,  the  equipment  vulnerability,  and  the  density  of 
facilities.  For  any  given  county  and  equipment  class,  the  expected  number 
of  equipment  failures  per  accident  is  proportional  to  the  amount  of 
carbon  fibers  released  and  the  density  of  facilities,  and  is  inversely 
proportional  to  the  mean  exposure  to  failure  for  the  equipment.  The 
actual  computation  of  X  is  done  by  summing  up  contributions  from  each 
county  in  the  U.S.  and  from  each  equipment  class.  The  mechanics  of 
these  computations  and  the  determination  of  the  probability  distribution 
of  the  number  of  failures  are  presented  in  Appendix  A. 

The  second  set  of  relationships  links  the  mean  and  standard  devia¬ 
tion  of  the  dollar  loss  in  a  single  accident  to  the  parameters  of  the 
distribution  for  the  number  of  equipment  failures  in  an  accident.  These 
relationships  are  based  on  standard  formulae  for  conditional  expectation, 
and  they  can  be  found,  for  example,  in  Parzen,  E.,  Stochastic  Processes, 
p.  55.  The  equations  imply  that  the  expected  value  of  L,  the  total  dollar 
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TABLE  5-1 


GLOSSARY  OF  SYMBOLS 


Mean  outside  exposure  to  failure 

Number  of  equipment  failures  in  an  accident 

Total  dollar  loss  in  a  single  accident 

Dollar  loss  resulting  from  a  single  equipment  failure 

Expected  number  of  equipments  of  type  j  that  fail  given 
an  accident 

Probability  that  i  pieces  of  equipment  fail  in  an 
accident 

Total  dollar  loss  annually  for  all  accidents 
Number  of  accidental  failures  involving  CF  nationally 
Expected  value  of  N^ 

Expectation 

Dummy  variable  to  denote  number  of  events 
Dummy  variable  for  dollar  loss 
Variance 

Variable  X  given  dummy  value  n 

Dummy  variable  for  dollar  loss  per  accident 


loss  in  a  single  accident,  is  proportional  to  X,  the  expected  number  of 
equipment  failures  in  an  accident,  and  that  the  variance  of  L  has  two 
terms,  one  which  is  proportional  to  X  and  one  which  is  proportional  to 
the  variance  of  the  number  of  failures  per  accident. 

The  final  set  of  important  relationships  links  the  statistics  of 
the  total  annual  dollar  loss  for  all  accidents  to  the  statistics  of  the 
dollar  loss  in  a  single  accident.  These  results  are  based  on  the  same 
type  of  conditional  expectations  relationships  referred  to  above.  The 
expected  value  of  the  annual  dollar  loss  is  proportional  to  the  number 
of  accidents  per  year  and  the  expected  value  of  the  dollar  loss  per 
accident.  The  variance  of  the  dollar  loss  per  year  is  approximately 
proportional  to  the  variance  of  the  dollar  loss  per  accident  and  the 
expected  number  of  accidents  per  year. 

To  convert  the  statistics  of  annual  dollar  loss  into  a  distribution, 
some  standard  statistical  methods  are  used.  The  results  obtained  and  the 
outcome  of  a  sensitivity  analysis,  are  presented  in  the  remainder  of  the 
chapter. 

5.2  COMPUTATION  OF  LOSSES  PER  INCIDENT 

The  computation  of  the  dollar  losses  per  automobile  accident  is 
performed  in  two  separate  steps.  In  the  first  step,  a  probability 
distribution  of  the  number  of  failures  contingent  upon  a  single  accident 
is  calculated.  In  the  second  step,  the  statistics  of  the  dollar  losses 
(rather  than  the  number  of  failures)  are  computed. 

An  analytic  methodology  was  developed  to  compute  the  distribution 
of  the  number  of  failures  contingent  on  a  single  fire  accident.  The 
methodology  is  based  upon  the  fact  that  for  a  given  county  and  equipment 
class,  the  number  of  failures  is  approximately  Poisson  distributed. 

This  is  due  to  the  extremely  low  probability  of  equipment  failure  at 
the  levels  of  exposure  typically  computed  for  automobile  fires.  Because 
the  dominant  variation  in  economic  losses  is  due  to  the  Poisson  failure 
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process,  this  methodology  does  not  require  detailed  modelling  of  release 
conditions  or  accident  locations.  As  shown  in  Appendix  A,  the  expected 
number  of  failures  per  accident  is  directly  proportional  to  the  geographic 
density  of  equipment  and  the  amount  of  fibers  released  and  inversely 
proportional  to  the  equipment's  mean  failure  level,  E. 

Implementation  of  the  Poisson  methodology  required  tabulation  of 
data  for  approximately  3,000  counties  in  the  United  States,  81  equipment 
categories,  and  several  possible  release  amounts.  To  handle  these  data, 
a  computer  model  was  developed  and  used  to  determine  the  distribution  of 
failures  contingent  upon  a  single  fire  incident.  Figure  5-1  describes 
the  logical  flow  of  the  model  and  its  extrapolation  to  the  national 
level.  As  explained  in  Appendix  A,  the  model  tabulates  a  mixture  of  a 
large  number  of  Poisson  random  variables.  There  is  a  separate  random 
variable  for  each  combination  of  county,  equipment  category  and  amounts 
released.  The  model  adds  up  the  probabilities  of  any  number  of  failures 
given  each  of  these  possible  combinations  and  weighs  them  by  the  appro¬ 
priate  conditional  probability  of  that  scenario.  The  result  is  the 
probability  that,  given  an  accident  in  some  county,  a  given  number  of 
failures  will  occur.  This  distribution  is  presented  in  Table  5-2. 

The  next  step  in  the  analysis  was  to  develop  the  distribution  of 
dollar  loss  given  an  accident.  The  mean  and  variance  of  the  dollar 
losses  per  accident  depend  on  the  statistics  of  the  number  of  failures 
and  of  the  dollar  loss  failure.  For  example,  if  there  were  five 
equipment  failures,  then  the  expected  value  of  the  dollar  losses  in  the 
accident  would  be  five  times  the  expected  value  of  the  dollar  loss  per 
accident,  and  the  variance  would  be  five  times  the  variance  of  the  dollar 
loss  per  accident. 

Formally,  we  used  the  computer-generated  values  of  X.,  the  expected 

J 

number  of  equipment  of  type  j  that  failed  given  an  accident.  On  an 
aggregate  basis,  the  X.'s  represent  failure  rates  for  the  given  equip- 

vJ 
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TABLE  5-2 


PROBABILITY  DISTRIBUTION  OF  NUMBER 
OF  FAILURES  GIVEN  AN  ACCIDENT 


Number  of  Failures 
0 
1 
2 

3 

4 

5 

>  5 
Mean 

Standard  Deviaion 


Probability 

.98 

1.8  X  10"2 
2  X  10"3 
2  X  10“'^ 
1  X  lO"'" 
4  X  10"5 
4  X  10“6 
0.022 
0.17 
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merit  classes  and  the  conditional  probabilities  that  any  given  failure 
is  of  type  j.  Thus 

X . 

Prob( Equipment  Type  j  Fails [Some  Equipment  Fails)  = 

j 

Using  this  probability  function  together  with  the  economic  loss  estimate 
described  in  Chapter  4,  we  developed  a  distribution  of  the  dollar  loss 
per  failure,  X^.  We  then  used  the  following  equations  to  find  the  mean 
and  variance  of  L,  the  total  loss  per  accident. 

EL  =  (EX  )  (EN) 

(5-1) 

Var  L  (EN)Var  Xq  +  (EXo)2  Var  N 

The  expectation  equation  simply  states  that  the  expectation  of  total 
dollar  loss  in  an  accident  is  equal  to  the  number  of  failures  times  the 
dollar  loss  per  failure.  There  are  two  terms  in  the  variance  expression. 
The  first  term  represents  the  variability  due  to  the  dollar  loss  per 
failure  distribution,  while  the  second  term  represents  the  variability 
in  the  number  of  failures  per  accident.  The  variance  equation  is  not 
exact  due  to  the  correlation  between  the  dollar  loss  per  failure  and 
the  number  of  failures.  The  precise  form  of  the  computations  is  presented 
in  Appendix  E.  Using  those  expressions,  we  derived  estimates  for  the 
dollar  losses  in  an  accident,  as  presented  in  Table  5-3. 

Although  our  methodology  does  not  permit  us  to  determine  the  precise 
distribution  of  dollar  losses  per  accident,  we  developed  upper  bounds  for 
these  probabilities  based  on  a  standard  result  from  probability  theory. 
This  result,  which  is  known  as  the  Chebyshev  inequality,  was  used  to 
determine  upper  bounds  for  the  probability  distribution  of  dollar  losses 
per  accident  as  well  as  upper  bounds  for  the  distribution  of  the  dollar 
losses  annually.  The  Chebyshev  inequality  (see,  for  example.  Feller, 
Introduction  to  Probability  Theory  and  Its  Applications,  Vol .  II,  p.  151) 
states  that: 

Prob  (L  ^  EL  +  ta(L))  1  1/t^ 
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TABLE  5-3 


STATISTICS  OF  ECONOMIC  CONSEQUENCES  FOR  A  SINGLE  ACCIDENT  (1993) 


Variable 

Symbol 

Variable  Name 

Expected  Value 

Standard 

Deviation 

Number  of  equipment 
failures  per  incident 

0.022 

0.17 

Dollar  loss  per 
failure 

$131 

$754 

L 

Total  dollar  loss 
per  incident 

$2.88 

$114 
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Thus,  the  probability  that  the  risk  is  more  than  10  standard  deviations 

_2 

above  the  mean  is  less  than  or  equal  to  10  .  Utilizing  the  Chebyshev 

inequality,  we  developed  Table  5-4*which  presents  upper  bounds  for  risk 
val ues. 

5.3  DERIVATION  OF  NATIONAL  LOSS  STATISTICS 

The  next  step  in  the  analysis  was  to  compute  the  national  risk  pro¬ 
file,  which  requires  only  a  knowledge  of  the  mean  and  variance  of  dollar 
losses  per  accident.  To  derive  the  national  risk  profile,  a  two-step 
procedure  was  employed.  These  steps  consisted  of: 

•  Computation  of  the  mean  and  the  variance  of  the  national 
risk  profile,  and 

•  Estimation  of  a  probability  distribution  based  on 
statistical  results. 

To  compute  the  mean  and  variance  of  the  national  risk  profile,  the  follow¬ 
ing  conditional  expectation  equations  were  utilized: 

E(L)  =  (EM)  EL 

Var(L)  =  (EM)  Var  L  +  (Var  M)  (EL)^ 

where 

L  =  Dollar  loss  per  accident 
U  =  National  dollar  loss 

M  =  Number  of  accidental  fires  with  CF  nationally 
EM  =  Expected  value  of  M 
EL  =  Expected  value  of  L 

As  noted  in  Chapter  3,  there  are  88  fire  accidents  annually  involving 
general  aviation  aircraft  using  carbon  fibers.  Assuming  that  the  number 
of  accidents  per  year  M  is  a  Poisson  random  variable,  then  EM  =  88,  Var  M 
88,  and  hence,  EU  =  253  and  =  $1,067.  These  statistics  are  summarized 
in  Table  5-5.  We  again  derived  an  upper  bound  for  this  distribution  based 

*A  second  version  of  the  inequality  used  only  for  the  first  entry  in  Table 
5-4,  states  that  Prob(L  It  EL)  1.  1/^ 
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TABLE  5-4 


UPPER  BOUNDS  FOR  THE  PROBABILITY  DISTRIBUTION  OF 

DOLLAR  LOSS 

PER  ACCIDENT  (1993) 

Dollar  Loss 

Upper  Bound  for  Probability 
that  Loss  Exceeds  this  Value 
Given  that  an  Accident  Occurs 

$  288 

10“2 

11,403 

10"^* 

114,000 

10"® 

1  ,140,000 

10"® 
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TABLE  5-5 


STATISTICS  OF  ECONOMIC  CONSEQUENCES  FOR  ALL 
ACCIDENTS  NATIONALLY  (1993) 


Variable 

Symbol 

Variable  Name 

Expected 

Value 

Standard 

Deviatior 

L 

Total  dollar  loss 
per  incident 

$2.88 

$114 

M 

Number  of  incidents 
per  year  (Poisson 

Di stribution) 

88 

88 

L 

Total  annual  dollar 
loss 

$253 

$1,067 
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on  the  Chebyshev  inequality.  These  results  are  presented  in  Table  5-6. 
The  national  risk  profile  is  depicted  graphically  in  Figure  5-2,  incorp¬ 
orating  the  Chebyshev  bounds  for  losses  in  excess  of  $50,000. 

5.4  SENSITIVITY  ANALYSIS 

We  next  examined  the  sensitivity  of  the  national  risk  profile  to 
input  assumptions.  Some  of  these  sensitivities  could  be  hand  calculated 
without  any  additional  computer  runs.  The  reason  for  this  is  that  the 
number  of  failures  per  accident  is  a  Poisson  random  variable.  Hence  the 
expected  value  and  variance  for  the  number  of  failures  are  approximately 
X  and  from  Equations  (5-1),  the  expected  loss  per  accident  is: 

X  E  X, 

and  the  variance  of  loss  per  accident  is  approximately  equal  to 

X(E  +  var  X^) 

As  an  example  of  a  sensitivity  analysis  using  these  equations, 
suppose  that  the  CF  amounts  released  in  an  accident  decrease  by  a  factor 
of  10.  In  this  case  the  expected  numbers  of  failures  for  the  various 
equipment  classes  would  all  decrease  by  a  factor  of  10,  while  the 
conditional  probability  of  dollar  loss  given  a  single  failure  would 
remain  the  same.  As  a  result  we  can  make  the  following  calculations 
for  the  loss  statistics.  Note  that  the  expected  national  loss  has 
decreased  by  a  factor  of  10,  to  $25. 

X  =  0.0022 
EL  =  0.29 
=  36 

EL  =  25.3 
a|-  =  338 

The  Chebyshev  inequality  results  are  tabulated  in  Table  5-7. 

We  also  examined  the  sensitivity  for  a  highly  conservative  scenario 
(Table  5-8)  in  which  amounts  released  were  increased  by  a  factor  of  ten. 
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TABLE  5-6 


CHEBYSHEV  BOUNDS  FOR  NATIONAL  RISK  PROFILE 


Annual  National 
Dollar  Loss 


Upper  Bound  for 
Probability  that  Loss 
Exceeds  Value 


10,923 

10"2 

106,953 

lO”** 

1 ,067,253 

10"6 
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APPROXIMATE  UPPER  BOUND  ON  NATIONAL  RISK  PROFILE  FOR  GENERAL  AVIATION  ACCIDENTS  (1993) 
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10,000  100,000  1  Million  10  Million 

Dollar  Value  D 


TABLE  5-7 


CHEBYSHEV  UPPER  BOUNDS  FOR  SENSITIVITY  ANALYSIS 
WHERE  RELEASE  AMOUNTS  DECREASE  BY  A  FACTOR  OF  TO 


Annual  Dollar 

Loss  for  Nation  (1993) 

Upper  Bound  for  Proability 
that  Loss  Exceeds  Value 

3,405 

10"2 

33,825 

lO"'' 

338,000 

10-^ 
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TABLE  5-8 


CHEBYSHEV  UPPER  BOUNDS  FOR  HIGHLY  CONSERVATIVE  SCENARIO 


Annual  Dollar 

Loss  for  Nation  (1993) 

Upper  Bound  for  Probability 
that  Loss  Exceeds  Value 

106,953 

10"3 

1  ,067,253 

10"5  , 

10,670,253 

10-7 
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The  same  sensitivity  analysis  applies  if  the  equipment  E  values  decrease 
by  a  factor  of  ten.  Note  that  the  loss  probabilities  increase  by  a  factor 
of  10  each  for  these  cases.  Thus,  for  example,  the  probability  of  exceed¬ 
ing  a  million  dollar  annual  loss  increases  from  approximately  10“®  to  10“® 
In  general,  if  the  amount  released  increases  by  a  given  factor  or  if  F 
values  decrease  by  the  same  factor,  then  the  probabilities  will  increase 
by  that  factor.  In  spite  of  the  highly  conservative  assumption,  the 
changed  values  in  Table  5-8  still  represent  low  probabilities  of  substan¬ 
tial  losses. 

5,5  SUMMARY  OF  RESULTS 

The  first  step  in  the  risk  analysis  number  was  to  project  the  number 
of  equipment  failures,  given  that  an  accident  occurred  somewhere  in  the 
U.S.  and  released  some  quantity  of  carbon  fibers.  The  expected  number 
of  failures  per  release  incident  was  extremely  small,  resulting  in  an 
expected  dollar  loss  per  incident  of  only  $2.88,  with  a  standard  deviation 
of  $114.  The  probability  of  an  accident  resulting  in  losses  exceeding 
$11,400  was  estimated  to  be  at  most  one  in  ten  thousand.  Then  based  on 
an  estimated  88  general  aviation  fire  accidents  per  year,  which  could 
potentially  release  CF  by  1993,  it  was  found  that  the  expected  annual 
loss  to  the  nation  as  a  whole  was  $253,  with  a  standard  deviation  of 
$1,067.  The  probability  that  the  national  loss  will  exceed  $107,000  was 
estimated  to  be  at  most  one  in  ten  thousand. 

The  sensitivity  of  these  results  to  several  input  parameters  was 
explored.  The  key  parameter  affecting  the  national  risk  is  the  amount 
of  carbon  fiber  which  could  potentially  be  released  in  an  accident. 

For  example,  decreasing  the  CF  release  quantities  by  a  factor  of  10  was 
found  to  decrease  the  national  risk  by  about  a  factor  of  10,  to  $25. 
Conversely,  increasing  the  CF  released  by  a  factor  of  10  would  increase 
the  expected  national  risk  to  about  $2,530.  The  chances  of  the  national 
losses  exceeding  $1.07  million  were  estimated  at  one  in  one-hundred 
thousand  for  this  scenario.  Hence,  this  highly  conservative  scenario 
also  results  in  low  probabilities  of  substantial  losses. 
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6.  CONCLUSIONS 


6.1  NATIONAL  RISK 

The  results  of  the  risk  analysis  indicate  that  the  potential  risks 
of  economic  losses  due  to  CF  releases  from  general  aviation  accidents 
are  relatively  small.  The  expected  national  risk  was  estimated  to  be 
only  about  $253  per  year  for  1993,  with  the  average  loss  per  incident 
being  on  the  order  of  a  few  dollars.  Furthermore,  the  chances  of  sub¬ 
stantial  losses  are  not  significant.  For  example  the  probability  of 
exceeding  $11,000  loss  in  one  year  was  estimated  to  be  about  1/100. 
Although  the  possible  consequences  of  a  single  accident  can  vary  greatly, 
depending  upon  whether  equipment  failures  do  occur,  the  likelihood  of 
such  a  failure  is  only  0.022  per  incident. 

It  should  be  noted,  however,  that  the  risk  estimates  are  subject  to 
uncertainty  from  a  number  of  different  sources.  The  assumptions  or 
uncertainties  incorporated  into  the  anlaysis  are  discussed  below.  Even 
when  sensitivity  analyses  were  performed  to  test  the  effect  of  these 
assumptions,  the  risks  were  found  to  be  reasonably  low.  For  example, 
the  likelihood  of  exceeding  $1  million  due  to  CF  releases  from  general 
aviation  accidents  is  only  10"^  even  if  the  amounts  released  are 
increased  by  a  factor  of  10. 

6.2  SUMMARY  OF  UNCERTAINTIES 

The  uncertainties  in  the  national  risk  estimate  may  be  analyzed 
by  considering  the  different  data  inputs  incorporated  into  the  model. 

The  chief  areas  of  uncertainty  are  the  fraction  of  fibers  released 
and  the  vulnerability  levels  of  electronic  equipment.  However,  even  the 
most  conservative  scenarios  in  our  sensitivity  analyses  indicate  that 
the  overall  national  risk  is  low.  Some  of  the  major  areas  of  uncertainty 
are  discussed  below: 

f  Carbon  fiber  usage  --  The  projected  usage  could  conceivably 
vary  by  a  factor  of  2  or  3  in  terms  of  CF  weight  per  air- 
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craft.  However,  such  variations  are  taken  into  account 
in  the  sensitivity  analysis  by  varying  the  amount  of  CF 
released  given  an  accident. 

•  Number  of  fibers  by  weight  --  The  present  report  assumes 
that  there  are  5  x  10^  single  fibers  per  kilogram  of  CF 
available  for  release,  based  on  previous  NASA  estimates. 
Although  this  value  is  subject  to  uncertainty,  5  x  10® 
represents  a  conservative  estimate. 

•  Fraction  of  CF  released  --  Recent  NASA  test  results 
indicate  that  the  2.9%  figure  used  in  our  base  analysis 
for  total  destruction  accidents  is  extremely  conservative, 
and  that  it  is  possible  that  no  more  than  0.1%  of  single 
fibers  by  weight  would  be  released.  In  any  case,  the  sensi¬ 
tivity  analysis  in  which  the  amount  of  CF  released  increases 
by  10,  covers  the  extreme  case  for  this  parameter. 

•  Accident  probability  --  The  estimate  of  340  accidents 
per  year  is  based  on  nine  years  of  historical  data.  A 
98%  Poisson  upper  confidence  bound  would  increase  this 
to  only  377.  Furthermore,  although  general  aviation 
activity  is  increasing,  there  is  no  statistical  evidence 
of  an  increase  in  the  number  of  accidents  occurring 
annual ly. 

•  Equipment  vulnerability  --  The  estimated  mean  failure 
levels  could  vary  considerably,  but  this  possibility 
was  addressed  in  the  sensitivity  analysis  for  the  con¬ 
servative  scenario  described  in  Chapter  5.  The  expected 
annual  losses  in  this  case,  also  assuming  a  ten-fold 
increase  in  CF  release,  were  about  $2500  for  1993. 
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•  Economic  losses  --  The  estimates  of  losses  per  equipment 
failure  are  subject  to  variations  between  facilities  and 
regions,  but  this  will  contribute  negligibly  to  the  over¬ 
all  uncertainty. 

In  summary,  the  sensitivity  analysis  indicates  that  there  is  very 
little  chance  of  substantial  losses,  with  a  "best  estimate"  expected 
annual  loss  for  1993  of  $250.  This  level  of  risk  is  quite  small 
compared  to  the  direct  property  damage  resulting  annually  from  general 
aviation  accidents.  However,  it  should  be  noted  that  the  present 
risk  assessment  has  addressed  only  dollar  losses  due  to  equipment  failure 
in  the  civilian  sector,  and  does  not  quantify  other  categories  of  risk 
such  as  costs  of  protection  or  cleanup  of  equipment,  CF  releases  from 

I 

non-aviation  sources  such  as  incineration  of  sporting  goods,  possible 
environmental  damage  by  carbon  fibers,  or  impacts  upon  military 
operations . 
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APPENDIX  A 


METHODOLOGY  AND  SYSTEM  OF  EQUATIONS  FOR  GENERAL  AVIATION  RISK  MODEL 


A.1  INTRODUCTION 

This  appendix  presents  the  methodology  and  procedure  for  constructing 
the  risk  profiles.  The  methodology  applies  the  Poisson  process  to 
release  types  and  is  based  on  actual  calculations  of  probabilities 
rather  than  a  simulation.  Section  A. 2  presents  the  rationale  for  the 
methodology  and  Section  A, 3  the  procedure. 


A. 2  BACKGROUND 

There  were  several  characteristics  that  distinguished  the  general  avia¬ 
tion  analysis  from  the  air  carrier  analysis  also  performed  by  Arthur  D. 
Little.  First,  the  collection  of  detailed  locational  data  on  accident 
scenarios  (locations  of  accidents  relative  to  locations  of  facilities) 
was  not  feasible.  Second,  the  expected  number  of  failed  pieces  of 
equipment  per  release  was  extremely  small.  Nearly  all  were  substan¬ 
tially  less  than  one. 

Given  these  differences,  a  different  type  of  methodology  was  used.  The 
basis  of  the  methodology  is  the  computation  of  the  expected  number  of 
failures  given  a  release  for  a  particular  equipment  type.  The  equation 
for  this  is: 


E(dA|r) 

=  /  n(dA)  (l  -  e  ^ 


(I) 


A-1 


where 


r  =  Set  of  release  conditions 
Nq  =  Expected  number  of  failures  for  given  release 

dA  =  Increment  of  surface  area 

A  =  Surface  area 

n(dA)  =  Density  of  equipment  within  area  increment  dA 
E(dA|r)  =  Exposure  within  area  increment  dA  given  set  of  release  conditions  r 

E^(dA)  =  Mean  exposure  to  failure  for  equipment  in  given  area  dA 
(incorporating  transfer  functions) 

For  general  aviation  accidents,  the  amounts  are  very  small  and  E  tends 
to  be  a  great  deal  smaller  than  Eq.  For  example,  a  typical  exposure 
contour  for  a  general  aviation  release  showed  maximum  exposures  of  10^ 
f.s/m^,  while  most  E^  values  are  at  least  10'^. 

In  view  of  this  (1)  can  be  approximated  using  Taylor  series  as: 

Although  n(dA)  may  not  be  uniform,  we  can  compute  the  average  value  of 
Nq  (averaged  over  release  conditions)  for  a  given  release  amount  of 
carbon  fibers  by 

N  - /f(r)  /n(dA)  dr  =  /  f(r)  /  E(dA|r)dr 

o  J  J  E  (dA)  Sr  A  '  ' 

r  A  0 

where 

ITq  =  Expected  number  of  failures  averaged  over  all  releases 
r  =  Release  conditions 
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=  Reciprocal  average  exposure  to  failure  of  equipment  in 
°  the  county 

f{r)  =  Probability  function  for  release  conditions 

and 

^0/  f(r)  /  n(dA)  dr 

n  =  - A - -  (4) 

/  f(r)  j  E(dA|r)  dr 
r  A 

In  other  words,  h  represents  the  average  density  of  equipment  where  the 
averaging  is  over  locations  weighted  by  exposure  and  vulnerability  values  for 
the  range  of  possible  release  conditions  for  a  given  amount  released.  Because  of 
the  random  locations  of  accidents  and  random  directions  of  wind,  T\  can 
be  approximated  by  D,  the  average  density  of  equipment  in  the  county. 

If  it  could  be  demonstrated  that  the  largest  concentrations  of  fibers 
generally  occur  at  the  locations  of  densest  concentrations,  then  rT 
would  exceed  D. 

To  investigate  the  possibility  that  n>D,  we  looked  at  average  city 
population  densities  weighted  by  population  (i.e.,  the  density  of  the 

city  of  the  average  person)  and  average  county  population  densities 
weighted  by  population.  Since  the  numbers  were  comparable,  we  concluded 

that  average  city  density  is  approximately  equal  to  average  county  density. 

By  extending  this  relation,  we  assumed  that  the  density  at  any  accident 
location  could  be  approximated  by  the  city  density  and  hence  the  county  density. 
Therefore,  as  a  first  approximation,  we  assumed 

D:Tr  (5) 

We  also  note  that  for  a  given  release  amount,  by  integrating  over  the  area  of 
exposure 

/  E(dAlr)  dA  =  S 
A 

where  S  denotes  the  surface  integral  of  exposure  and  is  a  constant. 


That  is,  the  surface  integral  of  exposure  is  simply  the  number  of  fibers 
released  times  the  settling  velocity.  Hencej  no  matter  what  the  weather 
conditions  are,  all  fibers  contribute  the  same  increments  to  the  surface 
integral.  Hence 

/f(r)  J  E(dAlr)  dr  =  /  S  dr  =  S  (6) 
r  A  r 


Combining  (3),  (5)  and  (6)  the  average  number  of  failures  for  a  given 
amount  released  is: 


N  : 

o 


DS 

E 

o 


(7) 


This  equation  was  one  of  the  two  key  results  of  the  analysis.  The 
other  key  result  was  that  the  distribution  of  the  number  of  failures 
is  approximately  Poisson  with  mean  Although  this  would  follow 
directly  if  it  could  be  assumed  that  the  individual  failures  are  negligibly 
correlated,  this  assumption  did  not  appear  to  be  immediately  justifiable. 
The  following  represents  an  alternative  approach  to  a  justification 
for  the  Poisson  distribution. 

Consider  the  random  variable  for  the  number  of  failures  given  TI  . 

In  computing  Nq  there  were  two  types  of  averaging  performed.  The  first 

type  was  averaging  the  random  failures  given  a  release  (i.e.,  the 
r 

average  in  Nq).  The  second  type  was  averaging  over  release  conditions 
such  as  stability  class,  wind  direction,  etc. 
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Given  the  exponential  failure  law,  then  the  number  of  failures  given 
the  average  is  Poisson  with  mean  equal  to  nJ  and  standard  deviation 
equal  to  V  ,  The  total  variation,  oased  on  a  formula  for  conditional 

expectation  (See,  for  example,  Parzen,  Stochastic  Processes,  P.  55) 
is 


Var(No.  FailuresIN^)  =  E(Var  No.  Failures |Nq) 

+  Var(E  FailureslNo) 

=  ENJ  +  Var  Ho 

The  first  term  EN^  is  the  Poisson  variation.  The  second  term  is  the 
variation  due  to  release  condition  and  density  variations. 


We  performed  some  computations  to  assess  the  relative  influence  of 

each  type  of  variation.  Table  A-1  presents  examples  of  total  deviations 

for  various  values  of  EN^.  It  is  assumed  in  the  Table  that  Var  N^  is 

four  times  EN*",  that  is,  the  standard  deviation  due  to  release  conditions 
0 

and  density  variations  is  double  the  mean.  This  was  simply  an  arbitrary 
but  in  our  judgement  conservative  assumption  made  to  test  the  importance 
of  the  Poisson  variation.  Column  1  is  EN^,  Column  2  is  the  Poisson 
standard  deviation  which  is  equal  to  the  standard  deviation  of  the 
number  of  failures  assuming  a  Poisson  assumption  is  valid.  That  is, 
the  second  column  is  the  square  root  of  the  first  term  of  the  right 
hand  side  of  the  above  equation.  Column  3  is  the  square  root  of  the 
right  hand  side  of  the  above  equation. 
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The  table  shows  that  the  assumption  of  a  Poisson  Process  with  parameter 
Nq  has  virtually  the  same  variation  as  the  actual  process.  Most  of  the 
actual  expectations  were  substantially  below  the  values  in  the  table. 

The  highest  expectations  were  household  goods  for  the  New  York  City 
counties.  For  these  cases,  the  densities  were  on  the  order  of 
56,000  per  square  mile  and  the  Eg  value,  incorporating  average  threshold 
values  was  3.4  x  10^,  Thus,  the  maximum  ¥q  for  a  helicopter  total 
destruction  release  was 

N  =  56,000  mi”^  x  - — -  x  1-09  kg 

°  1609^ 


9 

X  10 


kg 


T  .032 


m 

sec 


3.4  X  10 


9 


f  *  sec 

2 

m 


=  .218 


No  other  category  except  telephone  exchanges  and  forklift  equipment  yields 
values  that  even  come  close, to  these  household  goods  values.  Further¬ 
more,  for  the  high  density  equipment  categories,'  the  densities  and 
hence  Ng  will  not  show  a  great  deal  of  variation  with  respect  to 
release  conditions. 

The  conclusion  of  this  analysis  is  that  the  process  can  be  approximated 
by  a  Poisson  process  with  parameter  Ng  as  determined  by  Equation  (7). 
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TABLE  A-1 


EXAMPLES  OF  VARIATION  OF  FAILURES 


EXPECTATION 

ENJ 

POISSON  STANDARD  DEVIATION 

TOTAL  DEVIATION 

.25 

.5 

.7 

.10 

.32 

.37 

.05 

.22 

,23 

.01 

.1 

.101 

.005 

.071 

.071 

.001 

.032 

.032 
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To  account  for  different  equipment  types,  it  was  noted  that  the  number 
of  total  failures  is  the  sum  of  individual  Poisson  processes. 


A. 3  POISSON  MODEL 


Let 


i  =  1 ,  . 
j  =  1 ,  . 
k  =  1,  . 


,  N  be  the  counties 

,  M  be  the  equipment  SIC  category  combinations 
,  R  be  the  cases  of  amount  released 


Let 


\'jk  " 


(8) 


where 

S|^  =  Surface  integral  of  exposure  for  release  type  k 
D^.j  =  Density  of  equipment  in  county  i 

Average  reciprocal  exposure  to  failure  for  equipment  j 
incorporating  transfer  function  (This  was  assumed  to  be  inde¬ 
pendent  of  county). 

The  is  the  parameter  of  the  Poisson  Process  for  equipment  type  j, 

county  i  and  release  type  k.  Then  for  all  equipment  types,  the  para¬ 
meters  for  the  Poisson  Process  (which  is  the  sum  of  individual  Poisson 

Processes  )  is 


Then,  for  general  aviation  CF  releases: 


P.  = 


Q.  = 


p(n).  = 


Prob(county  i)  ^  weighted  factor  of  local  and  itinerant  operations 
(see  Appendix  B) 

Prob (release  k) 


Prob(n  failures)  =  E  PjQi,® 

•  I  IK 

I  ,k 


■^k  (iikT 


n! 


(10) 


and  the  average  failure  rate  is 


X  =  E  P.Q,  A., 

.  ,  I  k  ik 
I  ,k 


(11) 


Because  all  of  the  ^  values  will  be  small,  the  calculations  of  the  pro¬ 
babilities  in  (10)  will  be  needed  only  for  a  limited  set  of  values.  In 
order  to  compute  conditional  risk  profiles,  probabilities  of  equipment 
types  given  a  release  needed  to  be  computed.  Bayes'  theorem  is  utilized 
for  this  computation  as  follows.  The  prior  probability  of  scenario  k  is 


If  n  failures  from  a  release  are  observed  then  the  posterior  probability 
of  scenario  k  is 

p(i,k|n)  -  - 

E  p(n|  i  ,k)  p(i  ,k) 
i,k 
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’’A 


Now  given  n  failures  from  one  release  under  scenario  i,  k,  the  pro¬ 
bability  of  any  one  being  type  j  is 


p(j  |n,  i  ,k) 


Thus,  given  n  failures  from  one  release,  the  probability  that  the 
scenario  is  i,  k,  and  the  failure  is  type  j  is 


P (j .  i  ,k|n) 


=  p(j  In,  I  ,k>  p(l  ,k|n)  = 


EP,Q,  e 


EP,Q^  e 
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Thus , 


p(j  |n)  *  E  p(j  ,  Ikin) 
i,k 


iyi\  °  *nk 

I  P.Q  e 

•  I  I  K  IK 

I  ,k 


(12) 


For  the  case  n  =  1,  it  is  seen  that 

p(j  h)  ~ 

X 


(13) 


where 


i,k 


P.Q,  X. 

I  k  ijk 


The  computation  of  the  risk  profile  is  based  on  the  expressions  for 
p(n)  and  pj-  The  p(n)  values  were  used  by  themselves  to  represent 
a  risk  profile  for  a  single  accident  and  were  used  to  compute  moments 
for  numbers  of  failures.  Due  to  equation  (13)  and  the  fact  that  multiple 

failures  were  unlikely  p.  was  used  to  compute  aggregated  cost  moments 

J 

rather  than  p(j In) . 
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APPENDIX  B:  ANALYSIS  OF  GENERAL  AVIATION  ACCIDENTS 


B.1  INTRODUCTION 

In  order  to  estimate  the  potential  risks  due  to  carbon  fiber  releases 
from  general  aviation  accidents,  it  is  necessary  to  quantify  the  proba¬ 
bility  of  an  accident  in  any  given  county.  This  chapter  presents  the 
model  and  the  associated  analysis  for  general  aviation  accidents.  The 
model  is  based  on  accident  and  operation  records  compiled  by  the  National 
Transportation  Safety  Board.  It  incorporates  two  probability  factors: 
the  probability  of  a  general  aviation  fire  accident  and  the  conditional 
probability  that  the  accident  will  be  a  given  type.  The  type  includes 
aircraft  class,  damage  category,  and  whether  the  accident  occurs  during 
cruise  or  on  or  near  the  airport.  Factors  such  as  weather  and  proximity 
to  airport  are  discussed  but  not  utilized  in  the  model.  This  appendix 
also  includes  a  discussion  of  how  the  accident  model  is  interfaced  with 
the  risk  profile  calculation.  The  profile  is  based  on  county  data  and 
thus  requires  the  conditional  probability  that  a  given  accident  occurs 
in  a  given  county. 

B.2  DATA  BASE  AND  IMPORTANT  VARIABLES 

The  analysis  of  general  aviation  accidents  was  based  on  National 
Transportation  Safety  Board  data  tapes  for  nine  years  between  1968  and 
1976.  The  tapes  include  accidents  of  U.S.  and  foreign-registered  civil 
aircraft  that  occurred  in  the  U.S.,  and  contain  over  400  items  of  infor¬ 
mation  for  each  accident  or  incident.  This  includes  detailed  information 
about  the  accident  conditions,  the  aircraft  involved,  the  weather,  the 
pilot,  and  the  airport  if  the  accident  occurred  in  the  proximity  of  one. 

In  analyzing  these  tapes  we  focused  upon  accidents  caused  by  fire  or 
explosion  and  accidents  in  which  fire  or  explosion  occurred  after  impact. 
3058  accidents  of  both  types  can  be  identified  for  the  nine  year  period, 
and  for  each  of  these  the  following  items  of  data  were  extracted  for 
further  consideration: 
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1 .  Fi le  number 

2.  Date  of  occurrence 

3.  Location 

4.  Type  of  aircraft 

5.  Number  of  engines 

6.  Type  of  power 

7.  Aircraft  damage 

8.  Fire  after  impact 

9.  Type  of  accident 

10.  Phase  of  operation 

11.  Altitude  of  occurrence 

12.  Conditions  of  light 

13.  Type  of  weather 

14.  Airport  proximity 

15.  Lateral  distance  from  runway  centerline 

16.  Terrain 

17.  Ceiling 

18.  Visibility 

19.  Precipitation 

20.  Location  of  fire 

21 .  Fi re  damage 


The  objective  in  analyzing  these  data  for  the  3,058  fire  accidents 
was  to  develop  an  expression  of  the  form  P^(X)  that  represents  the  probability 
of  an  accident  of  type  t  under  conditions  X.  In  examining  a  variable 
to  determine  its  impact  on  the  accident  probability  function,  two  criteria 
were  applied.  First,  is  the  accident  probability  significantly  influenced 
by  the  variable  in  a  manner  that  can  be  utilized  in  a  general  aviation 
accident  model?  Second,  does  the  variable  affect  carbon  fiber  dispersion? 


In  analyzing  the  data  for  the  fire  or  explosion  accidents  for  general 


aviation,  we  identifed  the  following 
probabilities  or  release  condition. 

1.  Type  of  aircraft 

2.  Phase  of  operation 

3.  Level  of  damage 


variables  that  can  affect  accident 

4.  Weather 

5.  Location  relative  to  airport 


The  types  of  aircraft  are  classified  by  number  of  engines,  type  of 
power  and  the  existence  of  fixed  or  non-fixed  wings.  The  number  of 
accidents  for  each  type  is  presented  in  Table  B-1.  In  analyzing  some 
of  the  other  accident  variables,  it  was  noted  that  aircraft  type  affects 
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TABLE  B-1 


GENERAL  AVIATION  ACCIDENT  TYPES 


INVOLVING  FIRE  (1968-1976) 


General  Type  of  Aircraft 

Fixed  Wing 

2795 

He  1 i copter 

24it 

G 1 ider 

1 

Ba 1  loon 

15 

Blimp 

0 

Dirigible 

0 

Rocket 

0 

Convert i plane 

0 

Gyroplane 

3 

Other 

0 

Number  of  Engines 

0 

15 

1 

2354 

2 

666 

3 

0 

4 

23 

Type  of  Power 

Reciprocating  Engine 

2912 

Turbojet  Engine 

32 

Turboprop  Engine 

51 

Turbofan 

3 

None (G 1 i der) 

15 

Turboshaft 

45 
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variables  such  as  extent  of  damage  and  phase  of  operation.  In  addition, 
the  different  types  of  general  aviation  aircraft  may  be  utilizing  differ¬ 
ent  amounts  of  carbon  fiber  composite.  The  model  consequently  utilized 
three  classes  of  general  aviation  aircraft  and  these  are  discussed  in 
the  next  section. 

The  phase  of  operation  affects  carbon  fiber  release  conditions  and 
extent  of  damage  three  ways.  First,  fires  may  be  more  intense  depending 
on  the  phase  of  operation.  Take-off  accidents,  for  example,  would  involve 
more  fuel.  Second,  cruise  accidents,  because  they  take  place  off  the 
ground,  may  result  in  releases  with  different  dispersion  characteristics. 
Third,  accidents  taking  place  on  or  near  airports  may  take  place  in  more 
densely  settled  areas  than  for  cruise  accidents.  This,  of  course,  greatly 
affects  the  economic  impact  of  a  release. 

Because  of  the  methodology  utilized  in  analyzing  the  risk  of  carbon 
fiber  usage  in  general  aviation  aircraft,  the  first  two  effects  above  of 
phase  of  operation  are  significant  only  to  the  extent  that  they  affected 
the  amount  of  carbon  fiber  released.  The  reason  for  this  is  that  the 
surface  integral  method  discussed  in  Section  1.2  requires  only  the  amount 
of  fibers  released.  The  third  effect  above,  due  to  density  of  facilities 
and  population,  is  an  important  part  of  the  model.  The  split  of  general 
aviation  accidents  between  cruise  and  on-airport  location  was  utilized 
to  determine  county  locations  of  accidents. 

The  four  phases  of  operation  identified  by  NTSB  include: 

•  Take  off 

•  Landing 

•  Cruise 

•  Static  or  Taxi 

The  actual  split  utilized  in  the  model  is: 

§  Cruise 

t  On  or  near  airport 

NTSB  recognizes  several  levels  of  damage  in  aircraft  accidents  and  inci¬ 
dents.  It  is  assumed  that  a  carbon  fiber  release  can  occur  only  if  the 
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level  of  damage  is  total  destruction  or  substantial  damage.  For  each 
level  of  damage  a  different  amount  of  carbon  fibers  released  was  assumed. 

Weather  is  an  important  variable  in  aircraft  safety.  In  classifying 
weather  conditions  for  accidents,  NTSB  recognizes  the  three  different 
classifications. 

•  VFR  -  Visual  Flight  Rules 

•  IFR  -  Instrument  Flight  Rules 

•  Below  Minimum 

This  classification  appears  to  be  a  sufficent  statistic  for  weather 
variables  in  terms  of  their  effect  on  accident  probabilities.  For  air 
carrier  accidents  it  was  noted  that  49%  of  the  relevant  accidents  occurred 
in  IFR  or  below  minimum  conditions.  The  probability  of  IFR  or  below 
minimum  weather  at  the  major  hub  airports  (weighted  by  operations)'  is 
only  11%  and  this  indicates  that  the  probability  of  accident  per  oper¬ 
ation  in  IFR  weather  is  substantially  larger  than  for  VFR  weather.  For 
general  aviation  accidents  and  incidents,  the  percentage  of  IFR  and  Below 
Minimum  accidents  is  a  great  deal  lower.  Of  the  3,058  accidents  and 
incidents  in  the  data  base,  only  16%  of  the  non-static  accidents  occurred 
in  IFR  or  below  minimum  weather.  This  may  be  due  to  the  lower  incidence 
of  general  aviation  flights  in  IFR  weather  but  in  any  case  the  percentage 
is  a  great  deal  closer  to  the  national  incidence  of  IFR  weather  than  the 
percentage  for  air  carrier  accidents.  We  concluded  that  the  weather 
variable  need  not  play  as  important  a  role  in  a  general  aviation  model. 

The  percentage  of  IFR  and  below  minimum  accidents  was  higher  for 
larger  general  aviation  aircraft  (greater  than  one  engine  or  non-recipro¬ 
cating  engine)  and  may  also  be  correlated  with  the  operational  phase 
(cruise  or  on  or  near  airport).  However,  the  effect  of  the  IFR  probability 
at  any  airport  is  limited.  For  air  carriers,  even  given  the  large  increase 
in  accident  probabilities  due  to  IFR  weather,  the  weather  factor  does  not 
create  substantial  adjustments  in  airport  accident  probabilities.  Because 
weather  does  not  have  a  substantial  impact  on  general  aviation  accident 
statistics  and  because  of  the  difficulty  involved  in  making  such  adjust¬ 
ments,  we  did  not  incorporate  weather  into  the  general  aviation  accident 
model . 
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The  final  variable  considered  was  accident  location  for  those  acci¬ 
dents  occurring  on  or  near  airports.  In  NTSB  reports,  distance  from  the 
edge  of  the  nearest  runway  and  lateral  distance  from  the  runway  were 
reported  for  accidents  taking  place  off  the  airport.  Distance  from  the 
edge  of  the  runway  was  classified  into  1/4  mile  intervals  up  to  1  mile, 

1  mile  intervals  up  to  5  miles,  and  a  single  category  for  beyond  5  miles. 
In  examining  the  distance  distributions  for  different  types  of  general 
aircraft  phases  and  planes,  there  are  three  significantly  different  dis¬ 
tance  distributions.  These  distributions  are  for: 

•  Accidents  during  the  take-off  phase 

•  Accidents  of  single-engine  planes  during  landing 

•  Accidents  of  other  than  single-engine  planes  during 
landing 

However,  the  risk  model  for  general  aviation  is  not  capable  of  utilizing 
the  fine  detail  associated  with  the  probabilistic  distance  distribution 
for  accidents  taking  place  near  the  airport.  For  this  reason  we  did  not 
utilize  distance  distributions  within  the  accident  model.  Accidents 
are  classified  as  taking  place  either  during  cruise  or  on  or  near  air¬ 
port. 


It  should  also  be  noted  that  because  of  the  large  incidence  of  cruise 
accidents  for  general  aviation,  a  third  location  variable,  altitude  was 
also  considered.  Table  B-2  presents  an  altitude  distribution  for  the 
accident  and  incident  for  which  altitude  was  recorded.  Although  there 
are  a  substantial  number  of  accidents  taking  place  in  the  air,  we  again 
did  not  consider  altitude  because  of  the  surface  integration  technique 
which  considers  only  total  amount  released  as  a  relevant  variable. 

B.3  ANALYSIS  OF  ACCIDENT  VARIABLES 

In  view  of  the  qualifications  on  location  and  weather  noted  in  the 
previous  section  the  only  variables  analyzed  for  the  purposes  of  model 
development  were: 

•  Type  of  aircraft 

•  Phase  of  operation 

•  Level  of  damage 
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TABLE  B-2 


ALTITUDE  OF  CRUISE  ACCIDENTS  AND  INCIDENTS  FOR 
GENERAL  AVIATION  AIRCRAFT  (1968-1976) 


Altitude  of  Occurrence 
(Feet) 

0  -  500 

500  -  1,000 

1,000  -  2,000 

2,000  -  5,000 

5,000  -  10,000 

10,000  -  20,000 

20,000  -  30,000 

30,000 

On  the  Ground 
Unknown 

Tota  1 


Number  of  Accidents 

58 

k] 

36 

75 

55 

20 

2 

1 

974 

455 

1717 
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In  addition  to  an  analysis  of  these  three  variables  in  terms  of 
accident  frequencies,  statistics  have  been  compared  for  the  beginning, 
the  middle  and  end  of  the  nine-year  period  to  assess  the  possible 
existence  of  a  time  trend. 

The  data  utilized  in  performing  the  analysis  are  presented  in  a 
series  of  tables  that  follow.  The  analysis  also  included  some  other 
data  which  were  omitted  from  this  report  because  the  variables 
analyzed  in  these  data  are  not  ultimately  utilized  in  the  model. 

These  variables  include  IFR-VFR  weather  frequencies  and  distance  from 
the  end  of  the  runway  for  various  phases  of  operation  and  aircraft 
type  categories.  Weather  data  are  presented  in  aggregate  for  each 
subset  of  the  overall  data  base,  and  the  distance  distribution  for 
accidents  taking  place  near  the  airport  are  presented  for  the 
overall  data  base.  These  data  are  presented  for  general  reference 
and  for  corroboration  of  some  of  the  conclusions  presented  in  the 
previous  section.  For  example,  it  is  noted  that  the  probability  of 
IFR  or  below  minimum  is  16%  but  is  somewhat  higher  for  multiple 
engine  and  non-reciprocating  craft.  This,  however,  can  be  explained 
by  the  high  incidence  of  on-airport  accidents  for  these  classifica¬ 
tions,  since  on-airport  accidents  show  a  much  higher  percentage  of 
IFR  or  below  minimum  accidents.  Weather,  however,  is  not  judged 
to  be  a  relevent  variable  for  the  purposes  of  analyzing  risk.  (It 
can  also  be  noted  that  accidents  taking  place  near  airports  comprise 
a  distribution  with  most  of  its  mass  at  small  distances.) 

Tables  presented  include  the  following: 

Table  B3  -  Aggregate  Data  Base  Statistics 

Table  B4  -  Distance  Statistics 

Table  B5  -  Statistics  for  Fixed  Wing  Craft 
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AGGREGATE  DATA 


TOTAL  PLANE 
. WITH. FIRE 


TABLE 

BASE  STATISTICS 

ACCIDENTS  OH 
INVOLVEMENT  . 


B-3 

GENERAL  AVIATION  (1968-1976) 

INCIDENTS  3050. 


LEVEL  OF  DAMAGE  .  . 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
M 1 NGH/N ONE/ UNKNOWN 


359. 

2680. 

19. 


PHASE  OF  OPERATION/WEATHER 

3T  AT  IC/UNKNOWN 
NONSTAT I C 


dEL  MI N/UNKNCWN/NOT  REPORTED 


48. 

3010. 


392. 
2527. 
91  . 


PHASE  OF  UPEHAT ICN/LEVEL 
TAKEOFF 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI NDR/NUNE/UNKNOWN 


LAND  I NG 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
.MINOR/NONE/UNKNOWN. 

CRUI SE 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI NOR/NONE/UNKNOWN 

STATIC/TAXI 

_  SUBSTANTIAL  DAMAGE., 
TOTAL  DESTRUCTION 
MI NQR/NONE/UNKNCWN 

UNKNOWN  _  . 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI NOR/NONE/UNKNOWN 


OF  DAMAGE 

622. 

90. 
53  1  . 

1. 

638. 

100. 

533. 

._5.,. 

1  717. 

140. 

-  1565. 

12. 

58. 

2a._ 

- -  -  -  -  29.  ' 

1. 

. . . . 23  .  _  __ 

22! 

0. 


REASON  FOR  FIRE 

ZA  U  S"6 
RESULT 


"19  IT 
2867. 
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TABLE  B-4 


DISTANCE  STATISTICS  FOR  GENERAL  AVIATION  FIRES 


On  Ai rport 
On  Seaplane  Base 
On  Hell  port 

On  Barge/Shlp/Platform 
In  Traffic  Pattern 
Within  1/4  mi le 
Within  1/2  mile 
Within  3/4  mi  le 
Within  1  mile 
Within  2  mi  1 es 
Within  3  miles 
Within  4  mi  1 es 
Within  5  miles 
Beyond  5  miles 
Unknown/Not  Reported 


Ta ke-off 

Landing 

209 

257 

1 

0 

7 

5 

0 

0 

148 

143 

83 

24 

38 

24 

13 

8 

25 

22 

16 

22 

2 

17 

4 

15 

2 

9 

58 

84 

16 

8 
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TABLE  B-5 


STATISTICS  FOR  FIXED  WING  CRAFT 


TOTAL  Pt-ANE  ACCIDENTS  Qf<  INCIDENTS  2795* 
»ith  pipe  INVCLVEMFNT 


LEVEL  OF  DAMAGE 

SUBSTANT I AL  CAN  AGE 
TOTAL  nCST»^UCTICN 
MI  NCf./NCNE/UNKNC  VSN 


PHASE  OF  OPERAT I CN/teE ATHER 

ST  AT  IC  /UNKNGV^N 
NONSTATIC 
IFF 
VFR 

BEL  MI  N/ONKNClikN/NOT  FEPCRTEl: 


phase  of  operaticn/le vel  cf  damage 


37  • 

2  758* 

385  . 
^287. 
86  . 


308  . 
2468  • 

IS  * 


T AKFCFF 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCT  ICN 
MI NCR/NONE/UNKNC WN 


564 

-  78^ 

505* 

1* 


LANDING 

SUeSTANT I AL  DAMAGE 
TOTAL  DESTRUCT  ICN 
Ml  NGR/NCNE/UNKNC  V^N 


-  58S* . . 

90. 

494. 

5. 


CRtI SE 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCTICN 
HJ  NLLR/NONE/UNKNCifcN 


116* 

1433. 

12. 


STATIC/TAXI 

SUBSTANTIAL  DAM ACE 
TOTAL  OESTRUCTICN 
MI  NCR/NCNE/UNKNC  t»N 


23. 

17. 

1  « 


UNKNOWN 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCTICN 
MI NOR/NONE/UNKNCteN 


20. 


1. 

IS. 

0. 


REASON  FOR  FIRE 

175. 
2620. 


CAUSE 

RESULT 
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TABLE  B-6 


STATISTICS  FOR  NON-FIXED  WING  CRAFT 


TUTAi-  PLANE 
HITE  PIPE 


ACCIDENTS  Ok  INCIDENTS 
INVGL  VE^'E^T 


263. 


LEVEL  OF  OA^'AGE 

SUQSTANTIAL  O^^AGL 
TOTAL  DcSTkUC r  ICN 
IWINCE/NCNE/UNKNCHN 


PEASE  OF  CPEkAT ICN/HEATHEH 

STAT IC /UNKNOWN 
NONE  TaT I C 


BEL  .MIN/UNKNCHN/NDT  PEPCRTEC 


1  1  • 

252. 

7. 

240. 

5. 


PEASE  OF  OPEPATICN/LEVE 
TAKEOFF 

SOBSTANTIAl  CAWAGE 
TOTAL  DESTRUCTICN 
MI NOk/NONE/UNKNCHN 

LANDING 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINCR/NONE/UNKNCHN 

CROI 56 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
mingr/ncne/unknchn 

STATIC/TAXI 

substantial  damage 

TOTAL  DESTRUCTION 
Ml NCfi/NONE/UNKNCHN 

UNKNOWN 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MI NCR/NONE/UNKNCHN 


LE  CAM  AGE 

3d. 

-  12* 

2C. 

0. 

4S>.  - - 

10. 

39. 

0. 

156  . 

24. 
132. 
_  .0* 

17. 

5. 

.  _  _ _ _ —12* 

0. 

3. 

. -  _ 0* 

3  . 

0. 


REASON  FOP  FIRE 

CAUSE 
RESULT 


16* 
247  . 
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TABLE  B-7 


STATISTICS  FOR  ONE-ENGINE  CRAFT 


TOTAL  PLANE 
mTh  FIRE 


ACCIDENTS  OR 
INVGLVERENT 


INCIDENTS  2354. 


LEVEL  OF  DAMAGE 

SUBSTANTIAL  DAMAGE  24S. 
TOTAL  OESTRUCTICN  2105. 
MINCR/NONE/UNKNCVN  0. 


PHASE  CF  OPERATION/teEATHER 


STATIC/UNKNOHN 

36  . 

NONSTATIC 

2316. 

IFR 

237. 

VFR 

202b. 

OEL  MIN/UNKNCAN/NCT  REPCRTEO 

53. 

PHASE  OF  OPERAT ICN/LE VEL  OF 

LAHAQE 

TAKECFF 

4  72. 

SUBSTANTIAL  DAMAGE 

69. 

TOTAL  DESTRUCTION 

402. 

MINOR/NONE/UNKNCWN 

0. 

LANDING 

-  - 415. 

SUBSTANTIAL  DAMAGE 

66  . 

TOTAL  OESTRUCTICN 

349. 

MI NOR/NONE/UNKNGAN 

G. 

CfxLISE 

14  12. 

SUBSTANTIAL  DAMAGE 

103. 

TOTAL  OESTRUCTICN 

1309. 

MI NCR/NONE/UNKNCUN 

- 

. 0. 

STATIC/TAXI 

34. 

SUeSlANTIAL  DAMAGE 

10. 

TOTAL  OESTRUCTICN 

24. 

MI NCR/NONE/UNKNC VN 

0. 

UNKNCtaN 

2  1. 

SUBSTANTIAL  DAMAGE 

1. 

TOTAL  OESTRUCTICN 

20. 

MI NOR/NONE/UNKNCAN 

0. 

REASON  FOR  FIRF 

CAUSE 

loe. 

RESULT 

2246. 
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TABLE  B-8 


STATISTICS  FOR  OTHER  THAN  ONE-ENGINE  CRAFT 


total  plane 

*1  Th  FIRE 


ACCIDENTS  OR 
INVOLVEMENT 


INCIDENTS 


70A. 


LEVEL  OF  DAMAGE 

bOESTANTIAL  DAMAGE 
total  DESTRUCTICN 
MI  NCR/NGNE/UNKNCt»N 


PhA3E  OF  OPERA!  I  CNyVnEATHER 

STAT  IC/UNKNClfcN 
NGNST  AT  I  C 
IFR 
VFR 

BEL  M I N/UNKNCmN/NCT  REPORTED 


PHASE  OF  OPERAT ICN/LEVEL  OF  DAMAGE 
TAKEOFF 

SUESTANIIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MINCR/NONE/UNKNCy»N 


1  and  I  N<t 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MI  NCR/NCNE/UNKN’C  WN 

CRUISE 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MINCR/NONE/UNKNCI«N 

STATIC/TAXI 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MI  NCR/NONE/UNKNCtkN 

UNKNOV^N 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MINCR/NCNE/UNKNCMN 


REASCN  FOR  FIRE 

CAUSE  fi-* 

RESULT 
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;  tn  ut  :  no  GJ  j  rv)  hj  Caj  n  tn 

OMO  t-inOD  furis-vj  lai**^* 


TABLE  B-9 


STATISTICS  FOR  RECIPROCATING 

ENGINE  CRAFT 

TOTAL  PLANE  ACCIDENTS  Ort  INC 

IDENTS  2912. 

viiiE  fike  involvement 

■  ■■■  - - -  ■■ 

LEVEL  OF  DAMAGE 

SU.^STA^4T  I  AL  CAMAC-E 

33  1# 

TLIAL  DESTPLCTiCN 

2565  . 

Ml  NCF/NCNC/UNKNCV.N 

16« 

PHASE  OF  OPEHAT  ICN/^tEATHER 

STATIC/UNKNOkN 

NCNSTAT I C 

IFR 

VFR 

UEL  M I N/UNKNCVN/NOT  REPORTED 

44« 
2666  . 

36  9. 
2419. 
80. 

PHASE  OF  OPERA T I CN/LEVEL  OF  DAMAGE 

TAKEOFF 

SUOSTANTIAL  DAMAGE 

TOTAL  DESTROCTICN 

MI  NCfi/NONE/UNKNCi»N 

5S4* 

82. 

511. 

1. 

LANDING 

SUaSTANTIAL  DAMAGE 

TOTAL  DESTRUCTION 

MI  NCR/NDNE/UNKNCt»N 

_ _ «  0  1 

90. 
49  7. 

4. 

CRUI SE 

SUBSTANTIAL  DAMAGE 

TOTAL  DESTRUCTICN 
MINOB/NONE/UNKNCMN 

16  53* 

132. 

1510. 

-11. 

STATIC/TAXI 

SUBSTANTIAL  DAMAGE 

TOTAL  DESTRUCTION 

MI NOF/NONE/UNKNCUN 

52# 

26. 

,  6  . 

0. 

CNKNCInN 

SUBSTANTIAL  DAMAGE 

TOTAL  DESTRUCTION 

MI  NCfi/NUNE/UNKNCikN 

22. 

1  . 

2  1  . 

0. 

REASON  FOR  FIRE 

177* 

2735  . 


CAUSE 

RESULT 
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TABLE  B-10 

STATISTICS  FOR  NON-RECIPROCATING  ENGINE  CRAFT 

- 

lOTAL  PLANE  ACCIDENTS  OR  INCIDENTS 

146. 

- 

VilTH  FIKE  INVOLVEMENT 

LEVtL  OF  DAVAGC: 

suestantial  cavace 

28. 

TOTAL  DtSIRUCTICN 

115# 

yiNCA/NONE/UNKNChN 

3  • 

PHASE  OF  OPERATICN/AEAThER 

ST  AT IC/UNKNCteN 

4  • 

NONSTAT J  C 

142. 

IFk 

23. 

Vj-R 

1Q8. 

BEL  M I N/UNKNCWN/NCT  REPORTED 

1  1  • 

PHASE  OF  OPEHATICN/LEVEL  OF  CAMAGE 

T AKECFF 

28. 

SUBSTANTIAL  DAI^AQE 

a« 

TOTAL  OESTRUCTICN 

20. 

MINOR/NONE/UNKNCHpN 

0. 

LANDING 

47* 

.  ....  .  . 

SUBSTANTIAL  DAMAGE 

iO. 

TOTAL  OESTRUCTICN 

36  • 

MINCR/NONE/UNKNCVN 

1  « 

CRtI SF 

64« 

SUBSTANTIAL  DAMAGE 

6  • 

TOTAL  OESTRUCTICN 

56  • 

MINOR/NONE/UNKNCVN 

Xm 

STATIC/TAXI 

6. 

SUBSTANTIAL  DAMAGE 

2. 

TOTAL  DESTRUCTION 

.  — .  ... 

MINCR/NCNE/UNKNC VN 

1  . 

LNKNC  kN 

1  • 

SUESTANTIAL  DAMAGE 

.  ..  .  ... 

Q. 

TOTAL  OESTRUCTICN 

1  « 

- 

MI NOR/NONE/UNKNC VN 

0. 

REASON  FOR  FIRE 

- 

CAUSE 

14« 

RESULT 

132  « 
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TABLE  B-11 


STATISTICS  FOR  THE  PERIOD  1968-1970 


TOTAL  PLANE  ACCIDENTS  OR  INCIDENTS  1068* 
•  ITt-  FIRE  INVOLVEMENT 


LEVEL  OF  DAMAGE 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTHUCIICN 

MI ncf/none/onkncan 


113. 

948  • 

7* 


PHASE  OF  OPERAT ICN/AEATHER 

ST  AT  IC/UNKNO\*N 
NONSTATIC 


BEL  MI  N/JNKNCI»N/NOT  REPCRTHD 


"  20- 

1046. 


PHASE  OF  OPEHAT I  ON/LEVEL  OF  DAM  AGE 


TAKEOFF 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCTICN 
MI NCR/NONE/UNKNC  »N 

LANDING 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI NCH/NONE/UNKNCAN 

CRUI SE 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINOR/NONE/UNKNCVN 


233, 


231  • 


574- 


STATIC/TAXI 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI  NOR/NONE/UN KNCVN 

UNKNOWN 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MI  NCR/NONE/UNKNC  t»N 


REASON  FOR  FIRE 

CAUSE 

RESULT 


ei  • 

1007- 


117. 
913. 
I  8  . 


31. 

202. 

0. 


37. 
191  . 

3. 


35. 

535. 

.  .  .4. 


10. 

7. 

0. 


0. 

13. 

0. 
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TABLE  B-12 


STATISTICS  FOR  THE  PERIOD  1971-1973 


TOTAL  PLANE 
»ITt-  FIRE 


ACCIDENTS  OR 
INVCLVEMENT 


INCIDENTS 


910. 


LEVEL  CF  DAMAGE 

SUBSTANTIAL  DAMAGE  115. 
TOTAL  DESTRUCTION  769. 
MI NCF/NONE/UNKNCWN  O. 


PEASE  OF  OPEHAT  I  CN/KiEATHFfi 

STATIC/UNKNOKiN 

NONSTATIC 

IFR 

VFR 

BEL  MIN/UNKNOI»N/NOT  FEPCRTED 


1  1  . 

899. 


119. 

747. 

33. 


PHASE  OF  OPERATION/LEWEL 
TAKEOFF 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTHUCTICN 
MI NCB/NCNE/UNKNC AN 

LANDING 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINCR/NONE/UNKNCAN 

CRUI SE 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINOR/NONE/UNKNCkN 

STATIC/TAXI 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINGP/NONE/UNKNC  AN 

ONKNC AN 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCT ION 
MI NCR/NONE/UNKNC  AN 


OF  DAMAGE 

178  . 

29. 

14E. 

1. 

.192. 

33. 

158. 

1. 

521. 

45. 

472. 

4. 


REASON  FOR  FIRE 

63. 
847. 


CAUSE 

RESULT 
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QDtOO  0^00 


TABLE  B-13 


STATISTICS  FOR  THE  PERIOD  197^-1976 


TOTAL  PLANE  ACCICENTS  OW  INCIDENTS  1080. 
knITE  FINE  INVOLVEMENT 


LEVEL  CF  damage 


EU^fETANTIAL  DAMAGE  131. 
TOTAL  DuSTPUCTICN  943* 
MI  NL  h/NONE/LNKNC  VkN  6« 


phase  UF  OPEPATI CN/VkEATHER 

ST  AT  IC/UNKNOVlN  17* 

NONSTATIC  1063. 

IHk  156* 

VFf;  £c7« 

UEL  MIN/UNKNCVN/NOT  NEPCRTED  4C. 


PHASE  OF  OPEHAT ICN/LEVEL 
TAKECFF 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTPUCT  ICN 
MI  NCP/NCNE/UNKNCIaN 

LANDING  . 

SUaSTANTIAL  DAMAGE 
TOTAL  DESTRUCTION 
MINCP/NGNE/UNKNCmN 

CPUI SE 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCriCN 
MI  NCfi/NQNE/UNKNC'AN 

STATIC/TAXI 

SUBSTANTIAL  DAMAGE 
TOTAL  OESTRUCriCN 
MINCP/NCNE/UNKNCIfS 

LNKNQ\ftN 

SUBSTANTIAL  DAMAGE 
TOTAL  DESTRUCTICN 
MINCP/NaNE/UNKNCi\N 


CF  DAMAGE 

21  1  . 

. . . - . .  ..30* 

18  1. 

0. 

. .  — . . ...2. IS*-- - - 

30. 

184. 

1  . 

.  622 .  ' 

60. 

558. 

.  —,4. 

28. 

1C. 

. . 17. 

1  . 

4  . 

..  _ J. 

3. 

0. 


REASON  f-OP.  FI  (it 

67, 
1013. 


CAUSE 
RE  SELT 
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Table  B6  -  Statistics  for  Non-Fixed  Wing  Craft 
Table  B7  -  Statistics  for  One-engine  craft 
Table  B8  -  Statistics  for  Other  Than  One-engine  Craft 
Table  B9  -  Statistics  for  Reciprocating  Engine  Craft 
Table  BIO-  Statistics  for  Non-reciprocating  Craft 
Table  Bll-  Statistics  for  the  Period  1968-1970 

Table  B12-  Statistics  for  the  Period  1971-1973 

Table  B13-  Statistics  for  the  Period  1974-1976 

The  first  statistical  tests  examined  differences  in  the  statis¬ 

tics  for  the  three  3-year  periods  in  the  data  base.  These  included 
chi -squared  tests  for  the  following; 

•  Phase  (Take-off,  landing,  or  cruise)  versus  3-year  period 

•  Phase/weather  versus  3-year  period 

•  On  or  off  airport  versus  3-year  period 

•  Distance  from  airport  versus  3-year  period. 

There  was  no  significance  in  any  of  these  tests  with  one  minor 
exception.  For  the  phase  and  weather  data,  there  was  an  increase 
during  the  nine-year  period  in  below  minimum  or  unknown  accidents 
for  weather.  (This  may  be  due  to  reporting  practices.)  The  test 
based  on  data  with  a  deletion  of  the  unknown  or  below  minimum  cate¬ 
gory  shows  no  significance  in  differences  between  the  three  3-year 
periods. 
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Once  the  validity  of  the  full  nine  year  data  base  was  established, 
the  model  was  developed  in  a  straightforward  manner.  For  each  class  of 
aircraft  an  empirical  probability  was  estimated  giving  the  conditional 
probability  that  the  accident  took  place  either  on  or  near  the  airport 
or  during  cruise.  For  each  of  these  phases  the  conditional  probability 
of  substantial  damage  or  total  destruction  was  estimated.  One  class  of 
aircraft  with  unusual  characteristics  was  non-fixed  wing  or  non-powered 
craft.  These  craft  showed  a  different  statistics  from  fixed  wing  powered 
craft  and  in  addition  will  be  utilizing  different  amounts  of  carbon  fibers 
in  1973.  Most  of  these  aircraft  are,  of  course,  helicopters. 

In  examining  the  statistics  for  non-reciprocating  engine  planes  and 
planes  with  more  than  one  engine,  some  differences  from  one-engine  re¬ 
ciprocating  craft  were  noted.  There  was  a  higher  percentage  of  on-airport 
accidents  and  a  higher  percentage  of  substantial  damage  accidents  for 
these  craft.  Both  of  these  classes  of  aircraft  represent  larger  planes 
than  single  engine  reciprocating  engine  craft  and  are  logical  candidates 
for  consolidation  in  any  classification  scheme.  In  fact,  there  are  no 
significant  differences  in  phase  of  operation  or  in  damage  probabilities 
conditional  on  phase  between  the  two  classifications.  Although  there  is 
a  great  deal  of  overlap  among  the  two  classes,  that  is,  most  non-recipro¬ 
cating  engine  craft  are  also  aircraft  that  have  more  than  one  engine,  the 
similarity  of  statistics  shows  a  certain  pattern  for  larger  general  avia¬ 
tion  aircraft.  Because  of  this  the  final  breakdown  of  aircraft  was  estab¬ 
lished  as  follows: 

•  Non-Fixed  wing  or  non-powered  craft 

•  Single  reciprocating  engine  aircraft 

•  Multiple  or  non-reciprocating  engine  aircraft 
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The  statistics  for  these  classes  are  presented  in  Table  A14. 

The  conditional  probabilities  of  aircraft  type,  cruise,  on  airport  and 
damage  were  based  on  this  table.  As  noted,  the  major  differences 
between  the  statistics  of  the  larger  aircraft  and  the  smaller  single  re¬ 
ciprocating  engine  craft  are  the  larger  percentages  of  on  airport  acci¬ 
dents  and  substantial  damage  accidents.  To  some  extent,  these  differences 
are  correlated.  That  is,  larger  general  aviation  aircraft  have  higher 
percentages  of  substantial  damage  accidents  because  on  aircraft  (or  near 
airport)  accidents  have  higher  percentages  of  substantial  damage  accidents. 
In  fact,  the  conditional  damage  probabilities  given  the  major  phase  cate¬ 
gory  (on  or  near  airport  versus  cruise)  were  very  similar. 

B.4  GENERAL  AVIATION  ACCIDENT  RATES  AND  FINAL  MODEL 

Table  B15  presents  accident  rates  and  total  aircraft  hours  and  miles 
flown  for  the  period  1969-1978.  It  is  noted  that  the  number  of  fatal  and 
total  accidents  has  remained  approximately  constant  since  about  1972  while 
total  activity  has  increased.  The  model,  however,  does  not  require  a  per- 
mile,  per-hour,  or  per-operation  accident  rate.  Total  number  of  expected 
accidents  per  year  can  be  allocated  to  the  various  counties  according  to 
some  measure  of  activity.  Thus,  despite  the  decrease  in  accident  rate 
that  has  been  observed  over  the  past  several  years,  the  number  of  accidents 
has  not  substantially  varied,  and  we  assumed  that  the  expected  number  of 
general  aviation  fire  accidents  is  3,058  divided  by  9,  or  340  accidents  per 
year,  25%  (corresponding  to  the  percent  of  fleet  that  carry  CF)  or  85 
would  involve  CF  aircraft. 

It  should  be  noted  that  the  same  manufacturers  who  produce  general 
aviation  aircraft  also  produce  the  aircraft  that  are  utilized  by  air  taxi 
and  commuter  carriers.  We  did  not  perform  a  separate  study  of  air  taxi 
commuter  accidents  as  they  do  not  represent  a  significant  number  of  opera¬ 
tions.  In  1978,  for  example,  the  number  of  operations  for  each  aircraft 
carrier  category  at  F.A.A.  control  tower  airports  were: 
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Air  Carrier 
Air  Taxi 
G/A 

Military 

Total 


10,063,259 

3,773,484 

50,798,779 

2,537,912 

67,173,434 


The  1990  forecast,  however,  from  the  Wharton  long-term  industry  and 
economic  forecasting  model ^  predicts  8.4  million  air  taxi  and  commuter 
operations.  In  order  to  develop  an  estimate  for  the  impact  of  air  taxi 
and  commuter  operations  on  the  national  risk  profile,  we  adjusted  upward 
the  340  accident  per  year  to  account  for  air  taxi  and  commuter  operations. 
This  adjustment  factor  was  estimated  to  be  .02  which  is  consistent  with 
the  total  estimated -air  taxi  and  commuter  operations  and  general  aviation 
operations  in  1993  of  12.5  and  293.5  million,  respectively.  Thus  the 
total  number  of  estimated  general  aviation  and  air  taxi  and  commuter  fire 
accidents  was  354  per  year  of  which  88  would  involve  aircraft  carrying  CF. 
This  procedure,  of  course,  assumes  that  the  accident  rate  for  air  taxi  and 
commuter  operations  is  the  same  as  for  general  aviation  operations  and 
that  the  resulting  geographic  distribution  and  release  statistics  are 
similar. 


In  developing  the  final  conditional  probabilities  of  damage  and  type, 
the  figures  in  Table  B14  were  adjusted  to  account  for  different  growth 
rates  of  the  different  classes  of  aircraft.  As  noted  in  Table  2-1,  which 
was  used  in  project  future  fleet  sizes,  the  annual  growth  rates  have  been 
6.3%,  4.8%,  and  8.6%  for  multi-engine  craft,  single-engine  craft  and  rotor 
craft,  respectively.  We  would  expect  relative  accident  frequencies  for 
the  three  categories  in  Table  B14 to  grow  at  approximately  the  same  rate. 
(Thus,  for  example,  the  relative  share  of  helicopter  accidents  should 
increase.)  The  table  of  final  conditional  probabilities  reflects  these 
growth  rates  using  1973  as  the  base  year  for  the  Table  B14  data. 


^FAA  Aviation  Forecast,  Fiscal  years  1979  to  1990  U.S.  Department  of 
Transportation,  page  35. 
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STATISTICS  FOR  FINAL  CLASSIFICATIONS 
GENERAL  AVIATION  FIRES  (1968-1976) 
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The  final  part  of  the  model  requires  allocation  of  the  total  number 
of  accidents  to  the  various  counties  in  the  United  States.  A  separate 
allocation  procedure  was  developed  for  cruise  and  on  or  near  airport 
accidents.  For  on  or  near  airport  accidents  it  was  assumed  that  the 
probability  that  an  accident  takes  place  in  any  given  county  is  propor¬ 
tional  to  the  number  of  general  aviation  operations  in  that  county.  For 
the  905  airports  which  now  have  FAA  control  towers  or  are  candidates  for 
control  towers  we  obtained  the  number  of  operations  for  general  aviation 
as  well  as  air  taxi  and  commuter  for  the  year  1977,  and  we  obtained  fore¬ 
casts  for  the  years  1978  through  1990.^  We  then  used  the  forecast  to 
determine  a  growth  rate  to  predict  the  operations  at  these  airports  in 
1993.  These  airports  in  1977  included  an  estimated  30%  of  all  general 
aviation  operations  and  100%  of  all  air  taxi  and  commuter  operations. 

The  remainder  of  general  aviation  operations  for  1993  was  derived 
from  estimates  of  general  aviation  operations  by  state  for  1987.^  Again 
an  annual  growth  rate  was  used  to  estimate  state  operations  for  1993. 

From  the  total  operations  of  the  905  airports  and  the  total  state  opera¬ 
tions  we  determined  for  each  state  how  many  general  aviation  operations 
will  not  originate  or  terminate  at  one  of  the  905  airports.  These  un¬ 
allocated  operations  were  allocated  to  the  countries  within  the  state 
according  to  population.  They  were  then  added  to  the  totals  for  the  905 
airports  by  county  to  estimate  the  total  operations  by  county. 

While  accidents  that  take  place  on  or  near  an  airport  are  generally 
related  to  the  number  of  operations  at  the  airport,  accident  rates  for 
cruise  accidents,  on  the  other  hand,  might  be  more  appropriately  deter¬ 
mined  by  the  total  number  of  miles  or  hours.  An  estimate  of  total  hours 


^  Terminal  Area  Forecast  Fiscal  Years  1970  to  1990,  U.S.  Department  of 
Transportation,  Federal  Aviation  Administration,  1978. 

^  General  Aviation  Forecast  1975  to  1987  State,  Regional  and  National 
Operations,  prepared  for  the  U.S.  DOT  FAA  Office  of  Aviation  Policy, 
Final  Report,  April,  1976. 
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or  miles  by  county  is  thus  needed  to  allocate  cruise  accidents  to  the 
various  counties.  We  were  not,  however,  able  to  obtain  any  data  on 
general  aviation  mileage  wtihin  each  county.  By  the  same  type  of  proce¬ 
dures  noted  above,  we  were  able  to  obtain  forecasts  of  local  and  itinerant 
operations  for  each  county  in  the  country.  The  time  per  operation  of 
itinerant  and  local  flights  are  38  and  11  minutes,  respectively.  (Specif¬ 
ically,  for  example,  each  itinerant  operation  takes  on  average  76  minutes 
and  consists  on  average  of  5.4  operations.)  We  then  assumed  that  the  air 
mileage  of  each  itinerant  and  local  operation  is  approximately  the  same 
as  all  other  operations  of  the  same  type,  and  that  each  operation  takes 
place  completely  within  the  one  or  two  counties  where  the  take  offs  and 
landings  occur.  These  assumptions  are  obviously  not  true.  They  are  con¬ 
servative,  however,  in  the  sense  that  aircraft  mileage  will  be  assigned 
wholely  to  the  originating  and  terminating  counties  where  the  density  of 
vulnerable  facilities  is  generally  greater  than  any  other  county  of  the 
operations  path.  By  making  these  assumptions  the  probability  that  a 
given  cruise  accident  takes  place  within  a  given  county  is  proportional 
to  38  times  the  number  of  itinerant  operations  plus  eleven  times  the 
number  of  local  operations. 

Thus,  in  summary  the  general  aviation  accident  model  is  as  follows: 

1.  There  are  48  cruise  accidents  and  37  on  or  near  airport 
fire  accidents  per  year  (plus  an  additional  estimated  2 
and  1  for  air  taxi  and  commuter)  or  a  total  of  88. 

2.  For  each  cruise  accident 

38  Itin.  +  11  Loc. 

Prob(County  i)  =  — -  -  vo  T4-- - -  y - 

'  ^  38  Itin^  +  11  Loc^. 

where 

Itin^  =  Estimated  itinerant  operations  for  county  i 
Loc^  =  Estimated  local  operations  for  county  i 

and  for  each  on  or  near  airport  accident 
Prob(county  i)  =  ^tin^.  +  Loc^. _ 

Eintin.  +  Loc. 

1  1 

3.  The  conditional  probabilities  of  phase  and  aircraft  type,  based 
on  Table  B-14  and  adjusted  for  growth  rates,  are  as  presented 
in  Table  3-1. 
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EQUIPMENT  DATA 
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Facility  Category  -  Number  of  Pieces  $  Cost 

Category  Number  Equipment  E  Factor  Index  **  Per  Piece 


O 

00 


I/) 

TD 

O 

o 

CD 


OO  +J 


to 


cu 

{/) 


o 

IC 


00 


O  00 

LD  o 

CM  r- 


cr» 

00 


o  ro 

O  r— 

00  r- 

crt 


cn 


o  o 

CNJ  ^ 

o 


o 

oo 


o 

00 


o 

LO 

CM 


di 

o 

O)  03 

‘  c 

r-  S- 
-C  Z3 

S  Lu 


o 

-M 

03 

S- 

<u 

c 

cu 

CD 


O 
O  -r- 
•P  XJ 
O  03 

s: 


sz  x: 
•  (U  u 

to  (d 
c 

fd 


c. 

(U 


to 

c. 


C- 

O)  (U  I/) 

■M  C-  +j  r—  o  . 

^:3a)3i3<— .rd$-o 

5:-P  ^4->  CLQ.r~  C-PtO 

Q-ECSEr-.r-CC 
CUCTO'r-OOo3EOO 
Q-UJCJ  C-CDCD  E  S-OCJ 
>0  Q-  00  CU 

^  , —  , —  (L) - h-  O 

CU  O  r—  (U  r-  Cn 
r-tOfdCfdS^Xh-TD 
CU^  £*»-  E  03CQD:;  rd 
l—^CO  — lOO  iD_ooc: 


CU 


c. 

o 

-p 

03 

(U 

c 

<u 

CD  - 


C- 
o 

-p 
03 

cu 
c 
cu 

sCD  . 
,  CU 
S-  CD  S-  I 

0  4-0 
-P  03  -P 
O  r—  O 


03 


CM  CO 


to  r-- 


CO 


O  1 


CM 


CO 


o 

o 

00 


04 

r— 

o 

CO  CO 

o  o 

03 

CO  CO 

o  o 

r— 

O 

o  o 

03 

03 

03 

03 

O 

o 

+ 

o 

o 

+ 

^  LD 

o  o 

1  1 

in 
f  .3 

lO  IX) 

O  CD 

o 

o 

o  o 
o  o 

o 

06 

to 

o 

lO 

CD 

UJ 

o 

UJ 

o 

UJ  UJ 

LT)  CD 

LU 

UJ  UJ 

CM  <M- 

+ 

UJ 

o 

’f* 

UJ  UJ 

o  o 

\ 

UJ 

1 

UJ 

00 

7E- 

1 

UJ 

CM 

LO 

CO  »— 

CM 

CM  1— 

CO 

Kt 

CM 

CD 

O 

+ 

UJ 

CO 

.OE+06 

00  CO 

o  o 
+  4- 
UJ  LU 

CD  ^ 

.  1 E+09 

.8E+09 

.4E+09 

5E+06 

5E+07 

3E+07 

00 

o 

+ 

LU 

O 

00 

O 

4- 

LU 

O 

00 

O 

4- 

LU 

O 

00 

O 

4" 

LU 

LO 

CM  CM 

r— 

r—  CM 

r- 

•  • 

1—  r-. 

CM 

LO 

CD 

CM 

CD 

QZ 

* 

1—  CU 

in 

CD  CD 

.  jC 

O'— 

1 

s- 

O 

. — s  to  (U 

S-  o 

in 

0) 

_c 

r-  C2  > 

p  in 

c 

> 

P 

o 

1 —  rd 

c  C 

rd 

•r~ 

03  4.  C 

o  o 

4- 

cu 

S 

.  E  H- 

CD  CD 

h“ 

o 

to 

4. 

4. 

^  O 

o 

O 

CU 

cu 

•r~ 

•r— 

•r- 

P 

p 

X  T3 

"O 

T3 

4- 

4. 

CO  rd 

(d 

rd 

o 

o 

Q_  QZ 

DC 

CxZ 

to 

c>o 

cu 

u 


o 

Cl, 


cu 

C- 


s- 

cu 

p 

c: 

cu  cu 
o  o 

O) 

<4-  C 
O  •r-’ 
p 
P  4- 

t/1  o 
o  tn 

a. 


"O 

c 

cu 

Q. 

Ql 

<c 


cu 

P  in 
03  cu 
o  u 

C  •!“ 
23  TD 
4-  C 


C-1 


Equipment 

Facility  Category  -  Number  of  Pieces  $  Cost 

Category  Number  Equi pment  E  Factor  Index  Per  Piece 


o 

o 

o 

o 

O 

o  o 

o  o 

o 

o 

o 

o 

LO 

CO 

LD 

O 

LO  00 

o 

00  O  00 

o 

o 

o 

00 

CM 

CM 

00 

CM 

00 

00  r— 

eg 

r-v 

CO 

■* 

LO  I—  LD 
r- 


o  o  o 

cr» 

CT> 

cr» 

CT» 

cr» 

CTk 

cr> 

CO 

o  o 

o 

CM  CM  CM 

o 

o 

o 

o 

O 

o 

O 

o 

CM  O 

ro 

r—  O  r— 

CM 

00 

00 

CM 

00 

eg 

CM 

o 

o  o 

o 

o  o  o 

o 

o 

o 

O 

O 

o 

O 

o 

1  + 

1 

+  4-  + 

1 

1 

1 

! 

1 

1 

1 

4- 

UJ  Uul 

UJ 

LU  UJ  LU 

UJ 

LU 

UJ 

LU 

LU 

LU 

UJ 

LU 

O  LO 

o 

O  O  CO 

o 

CO 

ro 

O 

00 

O 

r—  O 

LO 

CM  r—  r— 

r— ’ 

oo 

CM 

00 

r*~ 

oo 

CM 

0^  CO 

00 

CO  CO  CTt 

CTi 

CJ> 

CO 

CO 

cr» 

cn 

o  o 

O 

o  o  o 

o 

O 

O 

o 

o 

o 

o 

o 

4-  4- 

4- 

4-4-4- 

4- 

4- 

■4- 

4- 

4“ 

4- 

4- 

4- 

LU  LU 

LU 

LU  LU  LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

LU 

00  LO 

CO 

cr>  "stj-  cr> 

o 

O 

OO 

OO 

00 

oo 

OO  CM 

CM 

oo  LO 

, _ 

CM 

LO 

00 

00 

00 

4-J 


• 

4J 

s- 

to 

, _ 

•r— 

r— 

C 

<U 

r— 

o 

1 - ' 

r— 

Q. 

• 

CO 

r— 

OJ 

O  XJ 

CD 

o 

s- 

tf— 

o 

£ 

CO 

r— 

o 

LO 

cu  OJ 

• 

U. 

CJ 

OJ 

S- 

4J 

u  El 

o 

1 

O 

c 

rd 

L- 

o 

cu 

JO 

1/) 

rd 

3 

r— 

4-> 

c 

03  p  ^ 

• 

o  ^ 

fd 

1 — 

C 

4J 

CO 

cu 

•  CL 

3 

UttC 

JZ 

S- 

JZ  S- 

C 

o 

CD 

+-> 

s- 

r— 

•1 — 

C 

c 

> 

crc>o 

OO 

4J 

o 

o 

03  CU 

T3 

o 

C-U 

CO  1 _ 

S-. 

CO 

.  fT3 

1— 

jd 

E 

o 

o 

•p— 

OJ 

Lf- 

3 

E  N 

o 

—  03 

n3 

>) 

E 

o 

CJ 

OJ 

^  r— 

r— 

• 

E  -r- 

CD  B 

Lu 

CO 

on 

OJ 

CO 

OJ 

(J 

Li-  03 

o 

4- 

-C 

-J 

! — 

rd  4-> 

i- 

• 

-P  CO 

.  o 

"■ — ^ 

r— 

^ > 

1— 

o 

on 

4J 

co 

o 

OJ 

CD 

CU 

S- 

3  N — 

O  T- 

• 

•r- 

•r— 

c: 

•  -r- 

+-> 

c 

-u> 

+-) 

onr- 

• 

-M 

4J 

CL 

“O 

-C 

X 

JD 

X 

l— 

-a 

03 

CO 

c 

03 

•r“ 

S- 

-p 

o  1— 

'1-3 

03 

rd 

t/> 

E  X 

rj  ro 

o 

CO 

O 

DQ 

QC 

rd 

S- 

fd  T- 

o 

:s 

o 

03 

S-  03 

c 

CD 

3 

C 

O  CO 

<  nc: 

CO 

D. 

21 

CU 

O  QC 

:=>  Q 

O  1— 

00 

U- 

CO 

Q-  Q_ 

l—H 

rc 

CU**-! 

O  Q. 

CO  CP>  O 
r—  r-  C\J 


I —  CvJ  CO  LO  <X) 

C\J  CNJ  CM  (NJ  CM  CM  ^ 


t — 

cn 

ro 

c 

c 

•r — 

•P“ 

S- 

E 

3 

■P 

CD 

O 

>> 

1- 

03 

rd 

03  4- 

2 

. 

3 

JZ 

cl: 

CD  C 

3 

• 

C  rd 

00 

cl: 

OJ  5! 

O 


C-2 


Equipment  categories  21-22  from  SIC  2824;  23  from  SIC  3714;  24-23  from 
SIC  2344 
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Facility  Category  -  Number  of  Pieces  $  Cost 

Category  Number  Equipment  E  Factor  Index  Per  Piece 
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APPENDIX  D 


Index 
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2 

3 

4 

5 

6 
7 

8-9 

10-11 

12-13 

14-15 

16-17 

18-19 

20-21 

22-23 

24-25 

26-27 

28-29 

30-31 

32-33 

34-35 

36-37 

38-39 

40-41 

42-43 

44-45 

46-47 

48-49 

50-51 


DEMOGRAPHIC  DATA  INDICES 

Demographic  Data  Category 

Dummy  Variable  =  1  Per  County 
Area 

Population 
Families 
Hospital  Beds 
Air  Carrier  Operations 
Number  of  Subway  Cars 
Facilities.  Empl oyees  SIC  Code  1900 

Facilities.  Employees  SIC  Code  3573 

Facilities.  Employees  SIC  Code  3650 

Facilities.  Employees  SIC  Code  3660 
Facilities.  Employees  SIC  Code  3670 
Facilities.  Employees  SIC  Code  3710 
Facilities.  Employees  SIC  Code  4011 
Facilities.  Employees  SIC  Code  4830 
Facilities.  Employees  SIC  Code  4890 
Facilities.  Employees  SIC  Code  5310 
Facil ities,  Empl oyees  SIC  Code  5410 
Facil ities.  Empl oyees  SIC  Code  5600 
Facilities,  Empl oyees  SIC  Code  5700 
Facilities,  Employees  SIC  Code  5900 
Facilities,  Employees  SIC  Code  6020 
Facil ities,  Empl oyees  SIC  Code  6100 
Facilities,  Employees  SIC  Code  6200 
Facilities,  Employees  SIC  Code  6300 
Facil i ties,  Employees  SIC  Code  7370 
Facilities,  Employees  SIC  Code  7391 
Facilities,  Employees  SIC  Code  8060 
Facilities,  Employees  SIC  Code  8220 
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APPENDIX  E 

DETAILS  ON  VARIANCE  OF  DOLLAR  LOSS  PER  ACCIDENT 


To  determine  the  dollar  loss  statistics  given  an  accident,  it  is  necessary 
to  condition  the  calculation  on  the  number  of  failures.  Thus  we  have 
the  following  conditional  expectation  formulae  for  the  first  and  second 
moments  of  dollar  loss  L  given  and  accident.  In  the  two  summations, 
the  variable  i  represents  the  number  of  failures  and  p(£)  represents 
the  probability  of  £  failures. 

EL  =  P(jl)  E(L|t) 
and 


EL2=  J,  p(il)  E(L2|t) 
and 


VAR  L  =  EL2  -  E(L)2 


where 


E(x|£)  represents  the  expected  value  of  variable  X  given  £  failures 
Because  expectations  are  additive,  we  have 
E(L|M  =  JlE(X^lii) 
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and 


EL  =  *P(t)  E(Xjl) 


=  EX  EN 
0 

where  X  represents  the  dollar  loss  per  failure  and  N  is  the  number 
of  failures  per  accident. 

For  the  variance  computation,  by  considering  the  individual  scenario 
probabilities,  one  can  derive  the  following  expression: 

EL2  -  p(j)  P,  (  iVar  (Xjl.k))  +  (xjl.k)^ 


where  now 

k  =  Amount  released 
i  =  County 

a  =  Number  of  failures 
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and  Q|^  are  the  scenario  probabilities  (See  Appendix  A)  and  the 

statistics  of  X  given  i  and  k  are  based  on  the  failure  rates  for  each 
0 

equipment  class  and  scenario.  An  alternate  expression  for  EL^  can  be 
obtained  by  considering  the  covariance  of  two  separate  losses  given  a 
failures.  The  expression  that  can  be  derived  in  this  case  is 

2  00  o 

EL  =  p(ii)  (£  var  (X^l^i)  +  E{x^\Si)  +  £(£  -  1)  cov^) 

where 


cov  =  covariance  of  two  separate  losses  given  a  failures 

The  approximate  expression  given  in  Chapter  5  (i.e.,  expression  5-1) 
follows  from  the  above  if  one  assumes  that  cov^  is  zero  and  that  the 
distribution  of  (X^U)  in  independent  of  i,  the  number  of  failures. 
These  assumptions  are  important  only  if  £  >  1.  Since  the  probability 
of  multiple  failures  is  very  low,  the  approximate  expression  in  Chapter 
5  is  virtually  identical  to  the  exact  expression  for  Var  L. 
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